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ABSTRACT 
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A decrease of one to six gauss (depending on magnetic 
field orientation) in the coupling constants of the CI^CC^ 
radical formed in X-irradiated zinc acetate dihydrate has 
been observed over a temoerature range of about -60°C to 
+30°C. Calculations of dipolar and Fermi contact interaction 
based on a model of internal rotation of the methylene grouo 
about the C-C bond have shown a small coupling constant de- 
crease on the order of 0.19 gauss over a temperature range 
of -150°C to +90°C. The major effect has been shown to be 
due to a spin relaxation mechanism. The effect was calculated 
using Monte Carlo techniques , the results of which were con- 
firmed by experimental data. 
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INTRODUCTION 



Study of the election spin resonance (esr) spectrum of 
single crystals damage by radiation have led researchers to 
the identification and structure of a number of paramagnetic 
species [Refs. 8, 11, 15, 13, 21, and 22]. Radicals mav be 
generated by irradiation of single crystals at nearly any 
temperature where the crystal remains intact. However, 
irradiation at one temperature will not necessarily produce 
the same radical or combination of radicals as that at an- 
other temperature. Hence, by selection of the appropriate 
temperature, it is possible to selectively study several 
radicals which may be generated in the same system. 

A considerably lesser amount of research has been done 
concerning the temperature dependence of the observed esr 
spectra for given paramagnetic species. Crawford [Ref. 4], 
in his study of the zinc acetate dihydrate crvstal and the 
paramagnetic species, CH 2 C0 2 , formed therein, noted a marked 
temperature effect in the sum of the couplina constants. 

The effect ranged from a one to six gauss change (for an av- 
erage coupling constant of about 45 gauss) over a temperature 
range of 298°K to 133°K depending on the external magnetic 
field orientation. Moss [Ref. 14] calculated the temperature 
effect in the esr spectrum of methyl radicals by attributing 
the effect to out-of-plane vibrations. However, he was unable 
to account for the magnitude of temperature coefficient bv 
applying his method to the benzene anion. In a study of the 
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hydrazine positive ion, N 2 H 4 + > Falle [Ref. 7] observed an in- 
crease in the ratio of nitrogen hyperfme splitting to oroton 
hyperfine splitting, j a n /a^ j , with increase in temperature. 

The effect was attributed primarily to a decrease in the cro- 
ton hyperfine splitting, a^ , as caused by torsional (rota- 
tional) displacement of an amino group about the N-N bond. 

This torsional motion causes an introduction of positive spin 
density at the protons leading to an increase in proton hy- 
perfine coupling, a^ , with increased torsional motion as ex- 
perienced with increasing temperature. Schrader [Ref. 17] 
has theoretically (based on the procedure of Schrader and 
Karplus [Ref. 18]) predicted the coupling constants for the 
methyl radical over a wide range of temperature from zero to 
500 degrees Kelvin, and has assigned temperature coefficients 
thusly which are in good agreement with Moss' results [Ref. 
14]. Zlochower, Miller, and Fraenkel (Ref. 23] also applied 
the theory of Schrader and Karplus [Ref. 18] to the methvl 
radical and observed experimentally a dependence of the pro- 
ton hyperfine splitting with temperature on the order of 
0.002 gauss/degree in a 50°C range extending above and below 
room temperature. Work done on the isotropic proton hyperfine 
interaction in the esr study of Tr-electron free radicals has 
shown that out-of-plane vibrations should have a contribution 
to the splitting as a result of direct overlap between the 
proton and the unpaired electron in the it orbital [Refs. 13 
and 14]. This contribution would be expected to increase with 
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increasing temperature leading to a corresponding decrease in 
the coupling constants with increasing temperature. 

Crawford [Ref. 4] gave no detailed exolanation of the 
marked shift in coupling constants observed for the CI^CC^ 
radical of zinc acetate. It becomes the dual ourpose of this 
work to verify the effect observed and then to develoo a model 
by which the phenomena may be explained. 
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II. BACKGROUND OF THE CH 2 C0 2 RADICAL 

Tolies, Crawford, and Valenti [Ref. 21] have recently 
published work on the eor (esr) study of zinc acetate dihy- 
drate single crystals. Almost simultaneously, Ohigashi and 
Kunta [Ref. 15] published a similar work. 

Both groups of authors have identified the CH 2 C0 2 radical 
to exist in X-irradiated single crystals of zinc acetate. 

The radical is stable for hours at room temperature with a 
half life of about ten minutes at 80°C. Figure 1 shows typ- 
ical esr spectra generated by the irradiated zinc acetate 
system. Coalescence of the typical four absorption line 
spectrum observed at low temperatures into the typical 1:2:1 
three line spectrum at high temperatures as indicated in Fig- 
ure 2 is evidence of internal rotation of the methylene group, 
CH 2 , about the C-C bond of the radical. This phenomena and 
determination of the barrier to internal rotation is discus- 
sed by Tolies et. al . [Ref. 21] and Ohigashi and Kunta 
(Ref. 15]. 

The crystal structure of zinc acetate dihydrate and the 
a, b, c cell axis system as related to an a*, b, c (x,y,z) 
orthogonal axis system within the crystal is shown in Figure 
3. The relationship between the two axis systems is qiven bv 
the transformations 
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Figure 1. Characteristic esr soectra at high and 
at low temperature. 
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Coalescence as 
Temperature 
Increases 



Figure 2. Coalescence of absorption lines of 
esr spectra with increase of 
temperature . 
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Figure 3. The zinc acetate dihydrate crystal and 
the orthogonal crystal axis system (xvz 
or a*bc) and its relationship to the 
crystalline axis system (abcj [Ref. 4]. 
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a = singx + cosBz 

b = y 

c = z 

where in Figure 3, x is also referred to as a*. 

The CI^CC^ radical is shown in Figure 4. The orientation 
shown is planar, each C-H group havmq assigned to it an 
orthogonal axis system , Y^ , Z^, or ^ 2 ' Z 2 w ^ erein the 

and Y 2 axes are coincident and fall along the p z orbital 
of the carbon (C 2 ) atom. These directions correspond to vec- 
tors u, v, and w in the a*, b, c crystal axis system. 

The following parameters of the radical as pre- 

sented by Tolies et. al . [Ref. 21] and Ohigashi and Kurita 
[Ref. 15] are pertinent to studies in this paper. 

Principal values of the hyperfine coupling tensors^ 

tensor -26.80 MHz 
-59.59 
-92.15 
H 2 tensor -26.41 
-58.98 
-92.78 



Principal values are designated to lie along the respec- 
tive C-H bond, the p orbital axis of the carbon atom in the 
n -- z 

CH 2 group, and in the CH 2 plane perpendicular to the respec- 
tive C-H^ bond as shown in Figure 4. 
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X 
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Figure 4 . 



Locations of the hyperfine axes in the 
CH 2 CC> 2 radical (planar configuration) 

[Ref. 4]. 
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Direction cosines with respect to the crystal axis 

. 1 
system 



H 1 


0 . 


. 547 


0 . 


,722 


0 . 


,423 




0 . 


. 768 


- 0 . 


,634 


0 . 


,088 




1 -0< 


. 332 


- 0 . 


,276 


0 . 


.902 


H 2 


0 . 


. 048 


- 0 . 


,081 


- 0 , 


.996 




0 . 


.768 


- 0 . 


.634 


0 , 


.088 




0 . 


. 638 


0 . 


. 769 


- 0 . 


.032 



The barrier potential to internal rotation of the 
methylene group about the C-C bond 

5.8 k cal/ mol 

The average rotational frequency factor 

12 -1 

4.4 x 10 sec . 



The direction cosines are expressed as row vectors. 
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III. PROBLEM FORMULATION 



To qualitatively understand the effect being studied, 
reference is made to Figure 1, page 15 which shows a typical 
first derivative esr spectrum of the CH 2 CO 2 radical at high 
temperature and at low temperature respectively. The esr ex- 
periment basically measures the interaction of the election 
magnetic moment with an applied magnetic field, H. Hvperfine 
structure is due to the free election interacting with the 
two methylene protons. The coupling constant, a H , is a 
property of the radical dependent on orientation in the mag- 
netic field and reflects the magnetic field of a given nu- 
cleus, , experienced by the electron. 

Figure 5 shows the splitting of energy levels in a maa- 
netic field due to electron and nuclei interactions for the 
CH 2 CO 2 radical. The four vertical lines represent allowed 
energy transitions giving rise to the basic four line spec- 
trum of Figure 1. 

When an oscillating field perpendicular to the external 
magnetic field H is applied, transitions are induced pro- 
vided the frequency, v, is such that the resonance condition 

hv = g0H 

is satisfied, where g is the spectroscopic splitting term and 
0 is the Bohr magneton. Hence it can be said 

AE. = hv. = g0H . 

1 1.1 

or 

H . “AE . 

1 1 
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Figure 5. Splitting of energy levels in a magnetic field due to electron 
and nuclei (H, and H«) interaction in the CH o C0 o radical. 
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Each resonance (absorption) line of the spectrum is a repre- 
sentation of a transition energy AE proportional to the 
magnetic field, H , causing resonance. Hence, the difference 
between two spectral lines m gauss is a representation of 
the difference in transition enerov which itself is a measure 
of the hyperfme interaction and is called the coupling con- 
stant. Experimentally the "distance" between absorption 
lines is directly measured in gauss to determine the relative 

coupling constants, a u and a , or the sum of the coupling 

H 1 H 2 

constants, a u + a^ seen in Figure 1. 

H 1 H 2 

As previously mentioned, Crawford [Ref. 4] in his work on 
zinc acetate dihydrate noted that at a given orientation of 
magnetic field, the esr spectrum showed a one to six gauss 
decrease in the sum of the coupling constants as the tempera- 
ture was increased from 133°K to room temperature (Figures 1 
and 2). In section I studies of temperature dependence of 
radical spectra have been discussed. These studies have con- 
tributed similar temperature effects to out-of-plane vibra- 
tions. Studies by Horsfield [Ref. 12] and Heller [Ref. 10] 
of methyl groups have shown internal rotations of the aroup 
to occur about a C-C bond. The temperature dependence of the 
CH 2 CO 2 radical has been attributed to motion of the CH 2 group 
consisting of rapid reorientation around the C-C bond above 
room temperatures being quenched at low temperature [Refs. 15 
and 21 ] . 
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It is the purpose of this work to develop a mathematical 
model, consistent: with experimental data, by which the temper- 
ature effect m the CH^CO^ radical may be predicted for anv 
given orientation of the zinc acetate crystal m a maanetic 
field. 

It will be seen m the discussion of results that the ob- 
served effect does not depend solely on one particular phenom- 
enon but rather is due to a combination of several phenomena. 
Hence the problem is formulated into four effects, rotation 
of the methylene group, hypercon^ ugation , flap or vibrations 
perpendicular to the plane of the radical , and nuclear soin 
relaxation . 

A. ROTATION OF THE METHYLENE CROUP 
1 . The Restricted Rotor Model 

The CH 0 (methylene group) acting as an internally 
restricted rotor representing the exchanging of the two pro- 
tons is the model used in explaining the temperature effect. 
The Hamiltonian for a restricted rotor is 

H = BL^ + ^(l-cos2p ' ) 

assuming the radical to be planar and the anqular momentum, 

L z , to be quantized along the Z axis (the C-C bond). In the 
Hamiltonian equation, B is a constant of the system, V is 
the activation energy or barrier to internal rotation, and 4> ' 
is the angle of rotation of the CH 2 rotor about the C-C bond. 

The Schroedinger wave equation for the restricted 
rotor may be solved by such standard procedures as given in 
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Eyrmg, Walter and Kimble [Ref. 6], page 358, to yie-ld wave 
functions which in the limit of the free rotor become: 



V = Ae 



im4> ' 



where A is a normalization constant egual to 



/2tt 



The set of basis functions are chosen to be 

, , im<f> ' / 1 , lmg ' 

|m-> = Ae Y = ( ) e 

/2tt 

where m runs from minus to plus infinity. 

Eigenvalues of the Hamiltonian yield the possible 
quantum mechanical energy levels of the system. These values 
are established by diagonalization of the Hamiltonian matrix, 
the diagonal values being the eigenvalues. 

The Hamiltonian matrix elements are determined 
through the following operations: 

On-diagonal matrix elements are represented as 



^m'|BL 2 +^|m> = <m'lH'im> = (Bm 2 +^)6 , 

1 z 2 1 11 2 mm ' 

where 6-0 unless m = m' . 

The diagonal elements are numerically evaluated from 
the knowledge of V, the barrier potential, and B, 
where 



V = 5.8 Kcal/mol 



h = Planck's constant 

I = the moment of inertia of the 
rotating group. 
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For convenience of calculation, 

B h 
h 87P 1 



where 



— 7 - = 5.05538 x 10 5 MHz amu A 2 . 

8 TT 2 

Off-diaqonal elements of the Hamiltonian matrix result 
from the perturbation term 

V i , 

— j COS<}> 

and are represented by 

<m' | - j cos4> ' | m> . 



Evaluation of these integrals shows that the value is 
zero unless 



m' = m+2 or m-2 . 

The non-zero off-diagonal elements are then expressed 
in terms of the basis function as 

im<}) ' 



-v r 

2 io 



2 t r -imij) 1 

cos2<j>' e 



/2tt 

Using the identity 



/2tT 



d<j>' 



0 , , 1 , i2(J) ' -i2(f) ' . 

cos2<j)' = j(e Y +e ) 



the integral is evaluated to yield 



8tt (2Tr<S (m-m' + 2) + 2tf6 (m-m’ -2 ) } 
or when m' = m+2 or m' = m-2 the off-diaqonal values 
become 

_V 

4 
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The calculated matrix elements are placed into the 



secular determinant which is diagonalized to yield the eigen- 
values and eigenvectors of the system. Solution of the deter- 
minant was performed by IBM 360 computer. 

2 . The Harmonic Oscillator Model 

The limits of a restricted rotor model are the free 
rotor at one extreme and the harmonic oscillator at the other. 
Considering the first several excited states above ground to 
be significantly populated at the temperatures used in this 
study, an approximately constant difference between levels of 
the restricted rotor was seen to exist (Section V-B) . 



gests the approximation that the system can be represented by 
a harmonic oscillator model, the energy levels of which are 
represented by 



where n runs from zero to infinity. If this approximation 
were to prove valid, further system analysis would be made 
easier . 



of E/h , hence some rearrangement of equation A2-1 is necessary 
to compare the two models. A solution of the Schroedinger 
equation for the harmonic oscillator by procedures of Eisberg 
[Ref. 5] shows 



The nearly constant difference in enerav levels suq- 




(A2-1) 



The restricted rotor model yields eigenvalues in units 




(A2-2 ) 



Equation A2-1 may then be rewritten 




(A2-3) 
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Direct substitution with proper units as used in section A 1 
yields 

E/h = 2 (n+y) ^18. 159208 (A2-4) 

from which eigenvalues may be determined m units of terahertz 
(THz) consistent with the units of the restricted rotor model. 
(One TKz equals 10^ Hz.) 

3 . Evaluation of Coupling Constants as a Function of 
Temperature 

The harmonic oscillator model was used to develop the 
values of coupling constants expected to be seen experimentally 
at any given temperature due to rotation of the CH 2 group. 

The effective spin Hamiltonian for the interaction 
of an electron and two nuclear spins in a maanetic field for 
the case of quenched angular momentum is 

d = BH *g * S - g 6 I *H + v S-T •! . (A3-1) 

^ r ^ n n 1 <f ~ 11 

1 = 1,2 1 = 1,2 

The first term represents the electron Zeeman energy, the 

interaction of an unpaired electron with the external maanetic 

* ■> 

field H. 6 is the Bohr magneton, S is the electronic spin 

equal to and g is the spectroscopic splitting tensor. The 

second term is the summation of the nuclear Zeeman energies 

where g is the nuclear g factor, 3 is the nuclear magneton, 

^n n 

t h 

and I is the spin of the 1 proton. The third term repre- 
sents the energy of hyperfine interaction between electron 

~ 7 ~ 

-> -> r 

spin, S, and nuclear spin, 1^. r iV is the hyperfine inter- 
action tensor. Equation A3-1 may be rewritten 
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(A3-2) 



» - eii-y-s - g n e n l ai 

1=1 ,2 n n 



where the parenthetical term is the effective maqnetic field 
seen by the i^ 1 proton. Solution of equation A3-2 leads to 
(among other parameters) the coupling constants of the spin 
system. 

When an electron and a proton are placed in the influ- 
ence of a magnetic field, the spins of each tend to line up 
parallel to the field creating dipolar coupling which is 
represented in the effective magnetic field term of equation 
A3-2. If the orientation of the magnetic field is chanced 
relative to the proton and electron as caused bv the rotation 
of the CH^ group, different dipolar coupling occurs, and the 
Hamiltonian solution yields different values of the couolina 
constants. Hence, the coupling constants are a function of 
the CH ^ angle of rotation, <f> ' . 

The Hamiltonian is a representation of the enerqv of 
the spin system. There will be a certain probabilitv P rf , , 
that a given energy level will be found at a given anqle of 
CH^ rotation, cf> ' . In addition, there will be a statistical 
distribution of electron spin population in the existing en- 
ergy levels as determined by Boltzman distribution factors. 

Based on the foregoing considerations the general re- 
lation to be solved to yield the average coupling constants 
for a given temperature is 



<A> 



T 




P , ,W d6 ' 

<p ' n Y 



(A3-3) 
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where <A>^, is the coupling constant averaaed over all orienta- 
tions for a given magnetic field orientation at a given tem- 
perature, T. , is the counting constant as a function of 

angle of rotation, <p ' , of the methylene groun about the C-C 

bond axis. P , , is the probability distribution factor, W 

tp n 

is the Boltzman statistical weighting factor. 

a. Evaluation of A,. 

<P 

An existing computer program, SSESR (Solid State 
Spectrum ESR) , was used as a starting point to calculate val- 
ues of A. ,. As originally written by Crawford [Ref. 4] , 

SSESR computes the esr spectrum of a radical trapped in the 
solid state assuming the angle c|> 1 to be zero, i.e. the planar 
model . 



The theory behind the SSESR program is based on 
the development of Box et. at [Ref. 2]. If the components 
of S perpendicular to H are ignored in equation A3-1 (a good 
approximation [Ref. 2]) the energy levels of the spin system 
can be written as 

E = 6gHS„ - g 6 El.. Ih.l (A3-4) 

^ H ^n n. ih' i 1 



where g is the observed g value, S fI is the electron spin com- 
ponent along H, and 1^ is the component of nuclear spin along 
hu , the vector along which the nuclear spin is quantized. The 

—y 

vector h^ may be written 




H 



(u. A 

l 



. a . +v . B . 3 • +w . C . 
li ill ii 



Y,) 



i q 3 
M n n 



(A3-5) 



In A3-5, A^ , B^ and are principal values of the hyperfine 

, —y ~y _ ~y , 

coupling tensors, T^. The vectors u^ , v^ , and w^ , are unit 
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vectors denoting the direction of the two sets of principal 
axes corresponding to A^ , B^ , and respectively, as seen in 
Figure 4. The terns ou , B^, and are the direction cosines 
of the external field, H, with respect to these axes. ^or 
the Ci^CC^ radical , 1 ^, and are each ±j. 

Incorporating all allowed combinations of spin 
states, solution of equation A3-4 yields 16 electron transi- 
tions in terms of 2AHgB n / c f B , represented bv 

U n n 



± l h K + )l 1 l h i(_)l 1 l h 2(+)l ± h 2 ( - ) 



(A3-6) 



where the 16 transitions are seen by the 2 choices of sign 

and where |h^ + j | and |h^_j I are the magnitudes of the vec- 

-* 11 
tors h. corresponding to S„ = + 75 - and respectively. The 

1 H Z Z 

computer program, SSESR, calculates the relative intensities 
and frequencies for the 16 transitions of the planar radical. 

Basic input data include the number of nuclei in 
the molecule or radical, maanetic field maanitude, direction 
cosines relating the nuclei principal axes to the crystal 
axis system, the principal values of the coupling tensors for 
each nucleus, the direction cosines of the maqnetic field, 
and the nuclear Zeeman frequency. Output yielded transition 
frequencies, transition intensities, and transitions in gauss 
from the center of the spectrum, all assuming an isotropic 
g tensor [Ref. 4] or a value of q = 2.0023, that for a free 
electron . 

Several modifications to the existing oroaram 
were made to incorporate the harmonic oscillator model. The 
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complete program, PREDICT II, appears in Appendix A. All 
computations were performed in double precision. 

(1) Subroutine HDCOS . The first modification 
involved an additional subroutine to compute the three direc- 
tion cosines of the magnetic field for a given input of 0 and 
<t> , the angles through which the crystal would be rotated in 
the magnetic field. This subroutine was called HDCOS and was 
based on the relations 



U 1 = 


-cos 


0 








c 

to 

II 


-sin 


e 


sin 


<l> 




U 3 = 


sin 


0 


cos 


4 > 


• 



(2) Subroutine AXROT . A second modification in- 
volved the subroutine AXROT, which computed the normalized 
direction cosines of the C-C bond axis using values of the 
direction cosines which relate the nuclei to the crvstal axis 
system as listed in section II. 

The vector, CH^ , representing that vector 
along the C-H, bond is expressed from values of the tensor 
(Section II) as 

x 1 = 0.547 
y ± = 0.722 
z^ = 0.423 

Likewise the vector CH 2 is represented by values of the H 2 
tensor 

x 2 = 0.048 

Y2 = -0.081 
z 2 = -0.996 . 
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Addition and normalization of these vector components yielded 
the direction cosines of the C-C bond with respect to the 
crystal axis system a*, b, c. 



rotates through angle <p ’ , the hyoerfine coupling tensors, 
and must change as a function of <p ' . Hence, a further 
modification by the inclusion of subroutine AXR was required. 
The rotation of any vector a number of degrees <J> * about a 
vector (a,b,c) is represented by the transformation matrix 
resulting from evaluation of the following matrix multiplica- 
tions : 



tor components of the C-C bond axis as computed in AX ROT , the 
resulting transformation matrix was applied to the and 
hyperfine tensors to yield the adjusted tensors as a function 
of <f> ' . These values then became input to the basic SSESR 
program. 



routine REFH was written to insure orthonormali tv of tensor 
values, the result of which became input to the AXR subroutine. 

With the modifications described it was pos- 
sible to calculate the stick spectra for any maanetic field 



(3) Subroutine AXR. As the methylene group 




where d = (l-c 2 ) 1/2 and a 2 + b 2 + c 2 = 1. 



Assigning a, b, and c the values of the vec- 



(4) Subroutine REFH. As a final refinement, sub- 
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orientation and any desired angle of rotation, <j> ' . From the 
stick spectra the coupling constants for any angle, <j> 1 , could 
be readily calculated by differences in the proper transition 
values . 



for the evaluation of the average coupling constants as a func- 
tion of temperature. 



solution of Schroedinger ' s equation for the harmonic oscil- 
lator as basically presented by Eisbera [Ref. 5] alona with 
a Lagrangian formulation was used to find the system's eiaen- 
functions , but the displacement variable was cj) ' rather than 
x as presented by Eisberg. 

The time independent Schroedinger wave equa- 
tion may be expressed as 



where V is the potential energy term, and I is the moment of 
inertia of the group. 



b. Evaluation of P , , and <A> m 

c p T 



Two additional modifications of SSESR were required 



(1) Angular Dependence of Wave Functions. The 




(A3-7) 



The potential energy may be expressed as 



V = |k<j> ' 2 . 



(A3- 8 ) 



The kinetic energy as 




(A3-9) 



Therefore, the Lagrangian (L) for the system 



is 



L = + 2 



(A3-10) 
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Solution of the Laqrangian yields the equation of motion 



4>* - 5<J>' = o 



(A3-11) 



which upon solution yields the frequency, v, as 



1 / K 

v = it/ t 

27T I 



( A3-12 ) 



To determine the actual value of the fre- 



quency, v, the following Hamiltonian was studied, 






(A3-13) 



where 



4> ' • 



21 



is derived directly from the classical mechanics 



expression for kinetic energy, 

KE = imv 2 

and L is angular momentum. 

Considering the angular momentum quantum 
mechanical operator, 

r - _*2 3 2 

L z " 11 3<J>' 2 



the Hamiltonian may be rewritten as 



n _ _h 2 3 2 , K 2 

21 3cf> ' 2 2^ 



(A3-14) 



which may be further rewritten using equation A3-12 as 

ti 2 3 2 



H = —. 



~i~2~ + V 2 2 7T 2 I 4) ' 2 



(A3-15) 



21 3 4 

Operating on the wave function, i p, with this 



Hamiltonian gives 

h 2 ^(cT) 



21 TV 2 " " + (v 27T iy = Elp (r) 



(A3-16) 



which is an identity with equation A3-14. Hence, rather than 
determine the frequency, v, by such means as, sav the IR 
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frequency between energy levels, it was possible to calculate 
K (and hence v) from equation A3-12 directly. 

The potential energy of the restricted rotor 

is expressed as 

V 2 

V = (1-cos 24>') (A3-17) 

where is the barrier potential observed experimentally as 

5.8 kcal/mol. Expanding the cosine term as a series vields 

„ V 2 V 2 , V 2 (24> ') 2 
v " 2 T~ 2 21 

+ higher terms. 

Neglecting higher terms , 

V = v 2 4>' 2 

which is just the potential energy of the harmonic oscillator, 
so that 

V = V 2 4>' 2 = 2-K4>' 2 (A3-18) 

and 



jK = V 2 = 5.8 kcal/mol. (A3-19) 

In dropping higher terms and approximating 
the harmonic oscillator the assumption has been made that 
significant values of energy will occur only with a few de- 
grees either side of 4> 1 = 0 and that as 4> ' approaches infin- 
ity, V approaches infinity as shown in Figure 6. 

Combining the results of equations A3-18 and 
A3-19, equation A3-7 is rewritten and transformed to give 



3 2 \p 



( 4 > ' ) 



34> : 



+ [ 



2 IE 



2I V A 



]<!> 



( 4 > ' ) 



= 0 



(A3-20) 
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Restricted Rotor 
Wave Function 




Harmonic 

Oscillator 

Potential 

Function 

and Wave 

Functions 



Figure 6. Comparison of restricted rotor and 
harmonic oscillator wave functions. 
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Defining the terms 



a = 


'2IV 2 1 


1/2 


( A3-21) 


Ti 2 




3 = 


2IE 

h 2 




(A3-22) 


and 


C = /a 4> ’ 


(A3-23) 



equation A3-20 becomes 
9 ( * i \ 

. 2 + = 0 

which with boundary conditions applied per Eisberg's proce- 
dures leads to the eigenfunction solutions 



iP U) = exp(-i; 2 /2)H nU) (A3-24) 

where H n ( ^ ) are the Hermite polynomials for the respective 
levels n = zero to n = infinity, n = zero representing ground 
state. It can be shown that these functions form an orthog- 
onal set over all space [Ref. 5], 

Combining equation A3-21 with equation A3-23 
and applying the proper units conversion factors , C was deter 
mined to be 



C = 3.77629 <j>' 
where <p ' is in radians. 



(A3-25) 



The probability, P , that a given state n 
will be found at an angle of CH 0 rotation, 4» ' , is 



P n(c) ( C) ^n ( 0 



where the complex conjugate of but since real 

wave functions are involved, 



P 



ur 




h 2 / \ 
n (?) 
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and the integral of equation A3-3 takes the form 



<A> t = 3.77629 



ft 00 _ 2 

\ A. ,e“ C H 2 , ,W dcj) ' 
♦) M <(> nU) n 



which when normalized, becomes 

90 



3.77629 

0 n , 1/2 
2 n! tt ' 



-90 



_ r 2 

e ^ H 2 . . W d<j) ' 
n (C) n 



(A3-26) 



Limits of integration were chosen based on the assumption 
previously mentioned that significant values occur onlv with- 
in a few degrees either side of $ ' =0. This assumption was 
justified by considering that for the ground state where 

H 0(C) = lf 




After about 0.27 radians or 16 degrees, P( z ) has decaved to 
— of its original value. After 64° it is expected that P, > 
will have negligible values, hence integration over the lim- 
its of -90° to +90° was considered suitablv close to those o f 
minus to plus infinity. 

(2) Temperature Dependence . As yet, no calcula- 
tion has been developed concerning temperature exolicitlv. 

The probability term derived above merely states the probabil- 
ity that a state n will be found at a criven angle of CI^ 
rotation <j> 1 . Boltzman statistics, however, when applied, 
show the statistical distribution of electrons between the 
possible energy levels, the greatest percentaae population 
occurring normally in the ground state level and decreasing 
exponentially with higher excited levels. 
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The population distribution in each energy 



level is expressed by the Boltzman weighting factor 



W = 
n 



g e 



-E /kT 
n' 



g e 



-E /kT 
ir 



Eg n e 

n 



-E /kT 
n / 



(A3-27) 



n 



where g^ is the degeneracy factor, is the difference be- 

t h 

tween zero energy and the n state energy, k is the Boltz- 
man constant and Z is the partition function. Since all 
levels in the model are doubly degenerate, (Section V-B) g 
cancels in the equation. 

( 3) Combining Coupling Constants, Angular and 
Temperature Dependence to Determine <A >,j 

Equation A3-3 may now be rewritten 
90 



<A> r 



f' 90 r 2 

= \ 3. 77629A , , T e ^ H 

t}_c ,n 4>' 



-E /kT 
n 7 



-90 



n=0 



n(0 



-dcj) ' (A3-28) 



Two subroutines were added to PREDICT II to 
perform the calculation. The first, AVGCC , evaluated the 
integral by numerical integration to yield the average cou- 
pling constants for each energy level, n, assuming 100% cop- 
ulation of each level. The second subroutine, AVBZM, applied 
the statistical factor to yield the average coupling constants 
as a function of temperature, <A> T » These two subroutines 
appear in Appendix A as an integral oart of PREDICT II. 



40 



B. HYPERCONJUGATION EFFECT ON COUPLING CONSTANTS 



Hydrogen atoms attached to 6 carbon atoms (carbon atoms 
adjacent to the carbon atom having unpaired spin) undergo 
interaction with the unpaired spin. This interaction is usu- 
ally termed hyperconjugation and is a maximum when the carbon- 
hydrogen bond and electron orbital are coplanar. Hydrogen 
atoms attached to an a carbon atom with an unpaired electron 
in a hybrid orbital also undergo interaction due to orbital 
overlap. This can be loosely considered a tvpe of hyper- 
conjugation and is the effect to be analyzed in this work. 

When the proton 2s orbital overlaps the carbon v z orbital 
there becomes a finite probability that the free electron 
will come in direct contact with the proton. This phenomenon 
is known as Fermi contact interaction and its effect is seen 
in the effective magnetic field term of equation A3-2. The 
Fermi contact interaction again changes the Hamiltonian val- 
ues and hence the coupling constants of the system. 

If the planar form of the C^CC^ radical is considered, 
a plane of symmetry exists and there will be no direct inter- 
action of the a carbon p z orbital and the hvdrogen 2s orbital 
as seen in Figure 7a. Spin-spin is the sole interaction in 
this configuration. At any other angle we expect direct or- 
bital interactions as seen in Figure 7b, the amount of Fermi 
contact interaction being determined by the electron density 
on the protons . 

The input coupling constants used in SSESR were those for 
the planar model. In order to adjust these values for Fermi 
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I 




I 

I 

Figure 7a. CH 2 orbital configuration in the planar 
form of the CH 2 C0 2 radical 




I 



Figure 7b. CH 2 orbital configuration after rotation 
of the CH 2 group through angle <j> 1 showing 
orbital overlap or mixing. 
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contact interaction (hyperconjugation) , a molecular orbital 
calculation [Refs. 4 and 19] using the CNDO approximation of 
Pople et.al [Ref. 16] was made to determine spin densities on 
the two protons for various angles of rotation, <J> ' . As seen 
in Section V, a plot of spin density versus angle of rotation 
yielded a hyperconjugation effect expressed by 

A . , = A. + B(0. 273xl0 - ^sin 2 2<t> 1 ) 

* ± 1 

where B = 1,420.4058 MHz, the hyperfine splitting constant of 
a hydrogen atom with a complete electron in the 2s orbital. 

This relationship was incorporated into the SSESR basic 
program as subroutine CCCA and used to adjust the computation 
of equation A3-3 for hyperconjugation as well as to determine 
the magnitude of the effect. CCCA appears as an integral 
part of PREDICT II in Appendix A. 

C. FLAP VIBRATIONAL MOTION 

An additional type of motion which can be termed as out- 
of-plane flapping motion of the CH 2 group may be considered 
as a possible source of contribution to the temperature de- 
pendence of the coupling constants, (figure 8). 

The effects of this type motion would be expected to act 
very similar to that discussed under the subject of hvper- 
conjugation except that in the flap motion, "bending" of the 
CH 2 group through an angle y is of interest. This action can 
be considered a further type of hyperconjugation effect lead- 
ing to mixing of the a carbon p z orbital with the hydrogen 2s 
orbitals but under a different mathematical relationship than 
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h4 







Figure 8 . Out-of-plane flapping motion of the 
CH2 group in the CH2CO2 radical. 
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that due to rotational motion. As with rotational motion, 
an increase of electron density on the protons (Fermi contact 
interaction) and a corresponding change in the coupling con- 
stants due to mixing of orbitals is expected for flap motion. 

Molecular orbital calculations based on the CNDO approx- 
imation (as used for rotational motion) were utilized to 
calculate the electron density for any angle of flap, y. 

From the results of the calculations, molecular energy val- 
ues were used to test the feasibility of this tvpe motion and 
electron spin density values to determine the functional de- 
pendence of the coupling constants upon the effect. 

D. SPIN RELAXATION EFFECT 

1 . General Relaxation Theory 

To observe electron spin resonance absorption there 
must be a population difference between the two spin levels 
between which transition occurs. Carrington and McLachlan 
[Ref. 3] discuss this criteria in rather easy to understand 
terms . 

The rate of absorption of energy from the radiation 
field (power) is 

H - Vae'W + Vbc'V V - nPAE (D1 - 1) 

where N is the population of enerav states E or E„. P „ 

a 6 a8 

is the probability of transition from state a to 6 (a|3) and 

is approximately equal to P. The difference in state popula- 
tion is n = N^ - Ng and follows a time dependent exponentiallv 
decaying function 

-2Pt 

n = n Q e (Dl-2) 
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Equations Dl-1 and Dl-2 show that resonance absorption dis- 
appears at equilibrium. 

Since N = N„ in the absence of a maanetic field, 
application of a field must inevitably require interactions 
between nuclei and their surroundings as thermal equilibrium 
is approached. This interaction is called spin-lattice re- 
laxation and results from changing spin orientation as mag- 
netic energy is transferred to other degrees of freedom. 

In addition, other processes are known to exist in 
spin systems due to nuclei-surroundings interaction which 
cause energy level shifts and hence an expected shift in cou- 
pling constants through solution of the effective spin Hamil- 
tonian of equation A3-1. These effects are called spin-spin 
relaxation, and through solutions of the Bloch equations 
yield correlation times representative of time the spin 
system remains in a given enercry configuration. 

2 . Application of Relaxation Theory 

Through a qualitative analysis of the physical pic- 
ture of the rotating methylene grouo, a relaxation effect was 
postulated. 

Consider as in Figure 9a the CH- methylene grouo 
oriented in a position having the maanetic moments of protons 
and H ^ precessing about effective maanetic fields reore- 
sented by vectors h + ^ and h + ^ consistent with the notation of 
Section A-3-a. When the two protons exchange sites as in 
Figure 9b, each will see a different effective magnetic field. 
There is then a finite probability that the protons will 
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Figure 9. Representation of the magnetic moments 
of nuclei and flipping from 

precession about one effective magnetic 
field to another due to rotation of the 
CH 2 group. 
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reorient or "flip" through an angle a into precession from 
about the original effective field to about the effective 
field now constituting their environment. The two limiting 
conditions occur when rotational frequency is very high corre- 
sponding to high temperature and when rotational frequencv is 
very low - essentially nonexistant - corresponding to quench- 
ing of motion at low temperature. Alternatively, it mav be 
said that at high temperature the frequency of rotation beincr 
fast will yield a high probability of flip (rapid flip) be- 
tween the two states. At low temperatures, the probability 
of making a flip will be small, and approach zero at the 
point where rotation can be considered to be cruenched. When 
the flip probability is low, significant relaxation effect is 
not expected. Hence a correspondingly greater average cou- 
pling constant is expected over that when flip probabilitv 
is high. 

It should be noted that Figure 9 is only intended to 
present a representative picture of the relaxation effect. 

In a real system, the probability of flip or reorientation of 
the magnetic moment must be considered at every increment of 
CH 2 rotation. 

3 . Calculation of Limiting Conditions 

In Section A-3-a Box's development for solution of the 
spin Hamiltonian was seen to yield transition frequencies and 
hence coupling constants for the CH 2 CO 2 radical. As derived 
from Box's relations the transition frequencies and hence 
coupling constant for low flip probabilitv (low temperature 
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limit) , are expressed by 



- h. 



- h. 



(D3-1 ) 



l l( + ) 1 1 “2 ( + ) 1 T 1 1 (-) 1 1 2 (-) 1 

where the results of both fields are assumed oriented in the 
same direction at quenched rotation. Solution of this rela- 
tion was obtained by a slight modification to Crawford's 
[Ref. 4] SSESR computer program. 



Coupling constants were obtained for the rapid flip 
(high temperature) limit through an approximation that at 
rapid flip conditions, the effective field experienced by 
each of the protons is an average of the two fields h^ and 
h^ experienced at quenched rotation. This approximation was 
expressed as 



l h l( + ) +h 2( + ) I + ' h l(-) +h 2(-) 



(D3-2) 



which was solved, again, by the modification of SSESR to 
yield the coupling constants for the rapid flip limit. 

The difference between the limiting values vielded a 
relative shift in coupling constants for the relaxation ef- 
fect at any given orientation of external maqnetic field. 

4 . Calculation of Intermediate Conditions 

Anderson [Ref. 1] has calculated frequency shift over 
a range of temperature but his calculations are based on an 
approximation that the motional Hamiltonian commutes with the 
time independent Hamiltonian, an approximation that does not 
apply to the model under consideration. To date a rigorous 
calculation without using Anderson's or other approximations 
has not been developed (Freed and Fraenkel [Ref. 9]). 
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As a result, a modified solution of the Bloch equa- 
tions in the presence of a random field usinq Monte Carlo 
techniques was used to calculate the frequency and hence, 
coupling constant shift for intermediate values of flip prob- 
ability (temperature) . 

A simple physical picture of the phenomenon calculated 
by the Monte Carlo technique is depicted in Figure 10 which 
represents a sphere with an orthogonal axis system X, Y, Z 
having its origin at the center of the sphere. At the origin 

— V- — V- 

is placed a proton. Vectors h^ and lie in the YZ plane at 
an angle a either side of the Z axis and represent the effec- 
tive magnetic fields h^ and h^ depicted in Figure 9. Assume, 
as an arbitrary starting point, that at the instant a maanetic 
field is turned on, the magnetic moment vector, M, of the pro- 
ton lies along the X axis with its tip falling on the surface 
of the sphere at polar coordinates (9 q, <(>q) . Let classical 
mechanics theory show that in the influence of a maqnetic 
field, the vector M will precess about h^ , tracing a oath 
along the sphere's surface from point one to point two. 
Finally, designate an increment of precession, A, such that 
after each increment a random decision as to whether or not 
to flip from precession about h^ to precession about f is 
made. The random decision is correlated to a selected prob- 
ability of no flip, P , which upon problem solution becomes 
correlated with a value of temperature. As time progresses, 
the magnetic moment will first flip from precession about 

— V - — >- 

to h 2 to h^ back and forth as shown after increment two, four, 
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Figure 10. Pictorial representation of the trace of 
the magnetic moment of a proton as it 
flips from precession about effective 

magnetic field h^ to , back and forth. 
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six, seven, eight, etc. The result is a random oath traced 
on the surface of the sphere. 

In the esr experiment an oscillating magnetic field 
along the X direction is observed. Hence, motion of the M^ 
component of the proton magnetic moment changing with time is 
the factor to be determined in the relaxation effect. If the 
probability of flip is zero, the M component will trace out 
a cosine-like curve. However, when a flip occurs and M 
changes its direction of precession, the M„ component deviates 
also from its original trace and begins tracing out a new 
cosine-like curve. As the M trace spreads over the surface 
of the sphere, it is expected that the M component will fol- 
low an approximate exponentially decaying cosine curve which 
can be represented by 

-a^t 

f,, , = e cosu a~t (D4-1) 

(t ) o 2 

where a^ is a decay constant, oJq is the frequencv with zero 
probability of flip, and a 2 is a frequency shift factor. 

Figure 11 is representative of the function, f( t p 
The Fourier transform of equation D4-1 (Figure lib) vields a 
single peak of frequency oj = characteristic of the flip 

probability and hence temperature. Plots of flip probability 
versus to may then be made for a qiven orientation of maanetic 
field from which coupling constant shift may be correlated 
with temperature. 

5 . Monte Carlo Calculation 

To calculate this effect, a random (Monte Carlo) tech- 
nique of generating the most probable function, f' (t) , 
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of magnetic moment. 



1 (oj) ( p ower) 




Figure lib. The Fourier transform of Figure 11a. 
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approximating that of Figure 11a was used. The comouter pro- 
gram was written by Tolies [Ref. 20] and appears in Appendix B. 

The magnetic moment M initially positioned at (9 q>4>q) 
was rotated initially about h^ then flipped from about to 
rotation about h^ , back and forth according to a first order 
Markoffian process. The problem is to find 



(w) 






0 



.. lU)t,, 

M , . . e dt 
x(t) 



where I ^ is a function of frequency and the integral is a 
Fourier integral giving the spectral density. 

Let A be an angle through which the magnetic moment, 
M, rotates about h^ prior to having the opportunity to flip 

—>■ — V t t t 

from state h^ to b .^ . The position of M after rotation by A 

-y 

about h 1 is expressed by the rotation matrix manioulations : 



'10 O’ 




/cosA -sinA 0 




'10 o’ 




/ sin0cos4> ' 


0 cosa sina 




sinA cosA 0 




0 cosa -sina 




sin0sin<J> 


i 0 -sina cosa 




O 

o 

I - 1 




( 0 sina cosa , 




cos 0 



cosA 

sinAcosa 



-sinAcosa sinAsina WsinGcostj) 

cosAcos 2 a+sin 2 a (-cosA+1) sinacosa sin0sin<f> 

COS0 

which yield a new point (0^,4>^) and hence a new value of the 



-sinAsina (-cosA+1) sinacosa cosAsin 2 a+cos 2 a 



M component. 

After a flip, rotation of A about yields another 
point (9 n / ( t' n ) expressed as above with each sina replaced by 
-sina. The rotational transformation is made after each in- 
crement of rotation, A. 
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To perform the Monte Carlo calculation a selected num- 
ber of points along the M trace separated by the value of A 

X 

* 

were chosen. For example, 36 points each separated by a A of 

-y -y 

10° would represent one cycle or one revolution of M about h, 
72 points would represent two cycles. A probabilitv, d, of 
having no spin flip after each increment A was chosen, a value 
of, for example d = 2 representing one chance out of two that 
the flip will not occur after increment A (rapid flipping) , a 
value of d = 100 representing 99 chances out of 100 that no 
flip will occur (slow flipping) . A random number chosen com- 
pletely arbitrarily was fed to a random number aenerator 
(RANDU) which produced in turn another completely random num- 
ber representing a condition of either "yes flip" or "no 
flip". RANDU was correlated to randomlv permit the system 
to flip in accordance with the desired flip probabilitv . . . 

2 or 100, etc. 

Starting with a point M Q ( 0 Q , 4> Q ) , rotation through 
A was calculated and a new value (0^, (f>^) was determined. At 
A the random number generated by RANDU indicated the system 
flipped or did not flip. If flip occurred, another value of 

■ y 

(0, <j>) was calculated based on rotation about h 2 . If no flip 
occurred, rotation was continued about h^ . After each incre- 
ment A, the decision was randomly made until the selected 
number of points was covered. Note that, especiallv for larae 
probability factors, such as 100, due to randomness of the 
generated yes/no decision, all points may be covered without 
any flip while on the other hand, several flips mav have 
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occurred. To obtain a most probable picture of the true ohys- 
ical phenomena, the calculation over the selected number of 
points was run over and over a selected number of times , the 
random number generator oroviding an ever increasing number 
of adjusted points along the M x (f' ( t ) ) curve. The final point 
values were then fitted to the damoed cosine function f ^ , 
equation D4-1, by a method of least squares. 

Output from the Monte Carlo calculation yields a trace 
of the damped cosine curve, the relative frequency, uiga 2 / the 
exponential decay constant a^ for the best fitting f^^, and 
a listing of deviations of the calculated f ' ^ curve from 
the best fitting f^ curve. 

6 . Probability as a Function of Temperature 

At the time, t, when the first flip decision is made 
(after A^) the probability, P , of no flip is 



1 



1 

d 



After 



2 TT 
A 



decisions, one cycle has elapsed and time 



t = — seconds. 

The probability, P^, of no flip has exponentially 



decayed to 



P = D 2,r/A 
n 



(D6-1) 



where 



D 5 




Hence equation D6-1 may be rewritten as 

-t/t 



P = e 
n 



(D6-2) 



where t is the time progressed and t c is a correlation time 
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representing the lifetime in a particular state before flio- 

ping to the alternate state. Using the relation 

_ InD 
D = e 

the probability of no flip, at t = ^ is 

P = exp(^ InD) = exp(-l/vT ) = exp(-l/t ). 
n a c o 



Hence , 



-A 

T c 2iTlnD 



(D6-4) 



where a basis of frequency of 1 Hz has been used. 

The results of Tolies et.al. [Ref. 21] and Ohigashi 
and Kurita [Ref. 15] show 



1 _ -V/RT 

= fn e 

T 0 



(D6-5) 



c 

where and V are the frequency factor and activation energy 
for the rotation of the CH 2 group. 

For given values of input parameters, A and d, T c 
may be calculated by equation D6-4. From these values, tem- 
peratures are determined by equation D6-5 resulting in a cor- 
relation between flip probability, temoerature , and frequency 
shift as calculated by the Monte Carlo computer program. 

7 . Fourier Representation of the Spectra 

The results of the Monte Carlo calculation yields 
valuable information about the predicted experimental spectrum 
representing the system. 

The symmetric Fourier integral of the curve, f , 
best fitting the calculated points of f' ^ is 



... • r 



• a. -a,t 

-10Jt 1 , 

e e cosa)Qa 2 tdt 
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Solution of this integral yields , 



-a. 



(t) 



a|+ (a>-a> 0 a 2 ) 2 



_ 1 _ 

a. 



1+ a^( w_a) 0 a 2 ] 2 



(D7-1 ) 



The Lorentzian line shape function, characteristic of damped 
oscillatory motion and characteristic of esr absorption spec- 
tra, is 

9 (u>) “ 1+T| (io-Wq) 2 ^ D7 

where T 2 represents relaxation effects. Comparing equations 
D7-1 and D7-2 the esr line shape parameters become 

Center frequency = u o a 2 

Half height = ^ = 2a. 

l 2 1 

Line width = = a. 

2 

These factors were used in analyses of the results of Sec- 
tion Y-F. 
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IV. EXPERIMENTAL PROCEDURES 



A. CRYSTAL GROWTH 

As described by Crawford [Ref. 4] single crystals of zinc 
acetate dihydrate were grown by slow evaporation of aqueous 
solutions. A mass ratio of about 2.5:1 water to zinc acetate 
dihydrate was used to form a nearly saturated solution, ali- 
quots of which were allowed to evaporate from watch glasses. 
Single crystals formed in 12 to 16 hours. 

B. CRYSTAL IRRADIATION 

Prior to irradiating a crystal, it was examined under a 
microscope to determine that the b axis (Figure 3) could be 
readily identified for proper orientation later in the esr 
spectrometer. A single crystal was placed on scotch tape and 
mounted in the open window of an X-ray diffractometer. Three 
such crystals could be irradiated simultaneously. Exposure 
was made typically for about two hours at room temperature 
using 50 KV, 25 ma , and a copper target. Immediately after 
irradiation, the crystals were placed in closed containers 
and refrigerated until use. The acetate radicals formed by 
the X-irradiation remained stable under refrigeration of zero 
degrees Centigrade up to one month. The signal may remain 
even longer, but no irradiated crystals were kept longer than 
this before study. 
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C. CRYSTAL ALIGNMENT IN THE MAGNETIC FIELD 



In order to vary the orientation of the magnetic field, 
a crystal was rotated through angles 0 and <f> as shown in 
Figure 12. The angles 0 and <J» are converted to the direction 
cosines of the magnetic field with respect to a*, b, c axis 
system through the transformations: 

Axis Direction Cosine 

a* -cos 0 

b -sin 0 sin <j> 

c sin 0 cos <(> 

Using these two angles it was then possible to describe 
any orientation in a three dimensional system. Since a 
slight change in the magnetic field orientation will produce 
significant change in the esr spectrum, precise orientation 
was necessary. 

A crystal to be studied in the esr spectrometer was 
placed on a quartz rod having both ends ground flat and par- 
allel to each other, one end for crystal mounting, the other 
for clamping to an extension rod connected above the resonance 
cavity to an angle indicator for determining the angle 0. 

The crystal was held on the quartz rod by a spot of stopcock 
grease acting as a glue. This procedure held the crvstal in 
place over the entire temperature range studied, -150°C to 
+100°C. 

Prior to placing the crystal into the spectrometer, the 
angle <p was established. As mentioned, careful orientation 
was critical. To establish <J> , a Bausch and Lomb microscope 
having a crosshair eyepiece over a calibrated revolving 
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Figure 12. Single crystal of zinc acetate dihydrate 
mounted on a quartz rod with the b axis 
coincident with the axis of the rod. 
Crystal is rotated by <j> and rod by 0 to 
establish magnetic field orientations. 
[Ref. 4 ] . 
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platform was used. The quartz rod was placed on the revolv- 
ing platform and one crosshair was aligned parallel to the 
long axis of the rod. The platform (and rod) was rotated 
through the desired angle <)> and the crystal was then aligned 
with the b axis parallel to the crosshair. The angle <p could 
be established to within ±0.5° by this procedure. 

The angle 0 was established after the crystal and rod 
were placed into the spectrometer resonance cavity. The end 
of the rod opposite to the mounted crystal was clamped to the 
extension rod leading to an angle indicator positioned di- 
rectly above the cavity. This indicator established angles 
to within ±1.0°. 

D. THE ESR SPECTROMETER AND TEMPERATURE CONTROL UNIT 

All data was taken using a Varian V-4502-13 esr spectrom- 
eter having a Fieldial V-Fr 2503 magnetic field regulator, a 
nine inch magnet, and a conventional reflection tvpe cavitv. 
Temperature control was attained using a Varian variable tem- 
perature control unit which had been calibrated for proper 
nitrogen flow to within ±1.5°C. 

During warmup and stabilization of the equipment, the 
angles <p and 0 were established. With both angles set, the 
spectrometer was fine tuned at a magnetic field of (3210 + 
40.00) gauss, the 40.00 being the initial setting of the 
Fieldial. The desired temperature was established using 
liquid nitrogen and nitrogen gas per operating characteristics 
of the temperature control unit. Ten minutes or longer was 
allowed for the crystal and cavity to reach temperature 
equilibrium. 
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Most data runs were taken at a set orientation of 9 and 
cf> , annotated as (^) , starting at a given temperature and vary- 
ing the temperature for each succeeding spectrum. Each time 
the temperature was changed, thermal equilibrium was allowed 
to be established prior to recording the spectrum. 

With all factors stabilized, first derivative spectra 
were recorded on an X-Y recorder over a 2.5 minute, 100 gauss 
sweep of magnetic field. Each absorption line of a spectrum 
was then determined and labeled using the field dial of the 
spectrometer. To avoid backlash error of the field dial, all 
lines were approached from the low magnetic field side when 
identifications were made. Line values were read to 0.01 
gauss. Measurements could be repeated to within ±0.03 gauss. 
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Y. RESULTS OF PROBLEM SOLUTION 



A. VERIFICATION OF SHIFT IN COUPLING CONSTANTS WITH 

TEMPERATURE 

The first objective in this investigation was to confirm 
the extraordinary shift in coupling constants with tempera- 
ture observed by Crawford [Ref. 4], In Section IV, the meth- 
od of recording experimental data was discussed. Generally, 
data was taken at a fixed orientation of magnetic field (|) 
while the temperature was initially set at -150°C and then 
raised by increments of ten or twenty degrees , spectra being 
recorded at each temperature increment up to about 90°C. 

Above 90°C the crystal began to dehydrate and the esr signal 
due to the C^CC^ radical rapidly broadened, decreased in 
intensity, and finally disappeared, no signal being observed 
at 100°C. This effect was probably due to molecular reorien- 
tation upon dehydration leading to destruction of the CH 2 CO 2 
radical . 

Figure 13 shows characteristic spectra obtained over the 
temperature range of -150° to +90°C. The particular spectra 
shown is that for a magnetic field orientation giving rise to 
two magnetically observable sites within the crystal. Two 
sets of coupling constants are obtained from spectra of this 
type, one for orientation (6), the other for orientation (_®) . 
When <J> is zero, the two magnetic sites are indistinguishable 
and spectra characteristic of that in Figure 1 are observed. 
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FIGURE 13 
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Data as recorded from the esr spectra is tabulated in 
Table I. From the difference in values (in gauss) of the 
absorption lines, experimental coupling constants were ex- 
tracted for a corresponding temperature. 

Figure 14 shows characteristic plots of the coupling con- 
stants versus temperature for given orientations of magnetic 
field, while the total shift in coupling constants between 
low and high temperature are tabulated in Table II. 

For all orientations of magnetic field, the total observ- 
ed shift in coupling constants was in close agreement with 
that observed by Crawford, verifying that the large shift 
does in fact occur for the CH 2 CC> 2 radical. Several addi- 
tional and very interesting phenomena were observed when the 
coupling constants were studied at increments of temperature 
as in Figure 14. Crawford's cursory observation of the tem- 
perature effect was based solely on the coupling constants 
determined at two points, + 25°C and -140°C. This investiga- 
tion confirms Crawford's observation but in addition shows 
that for all orientations of magnetic field, the temperature 
effect is confined to a narrow range of about +30°C to -60°C 
as seen by the characteristic "S" shaped plots. A decrease 
in the temperature coefficient was to be expected at low tem- 
peratures as the system approached quenched rotation of the 
CH 2 group, however the decrease at high temperatures above 
+30°C was totally unexpected and eventually led to formulation 
of the relaxation model described in Section III. Initially 
it was thought that the decrease in temperature coefficient 
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TABLE I. Experimental data. 



Field Temp 


Field 


Diff . 


Field Temp 


Field 


Diff. 


Orien. (°C) 


Dial ^ 


Taken 


Orien. (°C) 


Dial i 


Taken 




Reading 






Reading 




( 9 o - 150 


21.56 


39.81 


00 

+ 

'"N 

o O' 
<T\ 


21.18 


40.14 




31.91 






* 






51.12 






* 






61.37 






61.32 




-130 


21.55 


39.93 


+ 27 


21.12 


40.16 




31.95 






* 






51.19 






* 






61.48 






61.29 




-110 


21.47 


40.01 


+ 46 


21.08 


40.11 




31.88 






41.21 






51.13 






61.19 






61.48 




+ 60.5 


21.02 


40.10 


- 90 


21.41 


40.08 




41.13 






31.89 






61.12 






51.08 

61.49 




+ 81 


20.98 

41.01 


39.99 


- 70.5 


21.34 


40.14 




60.97 






31.95 

50.96 




+ 90.5 


20.94 


39 .95 




61.48 






40.96 

60.89 




- 50 


21.30 


40.18 


3 -so 

-70 








32.00 




2 

Footnote 


40.50 




50.88 






40.44 




61.48 




-63 




40.35 








-50 




40.25 


- 30.5 


21.25 


40.23 


-43 




40.15 




32.61 




-33 




39.85 




* 




-23 




39 . 31 




61.48 




-13 




38.30 








-11 




37.98 


- 11 


21.21 


40.19 


- 3 




37.65 




* 




+ 8 




37.20 




* 




+ 17 




36 . 75 




61.40 




+ 27 




36.59 








+ 51 




36.30 ' f 



Gauss field dial on spectrometer. 

2 . 

Field dial readings will not be listed for remaining data. 
3 

Data as corrected for Temperature recalibration. 

★ 

Unresolvable . 
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Field 


Temp. 


Differences 


Field 


TemD . 


Di f ferences 


.Orien . 




(°C) 


Taken 

j 


i Orien. 

t 


(°C) 


Taken 


j/120; 

!\ 0 J 


+ 60 . 5 


I 

36 . 39 


(180', 1 
Vo ! 


-90 


46.58 


+ 67 


36.25 


-50 


46.41 


i 


+ 75 


36 . 25 


\ 


-43 


46.20 




+ 90 


36.35 ! 




-33 


46.00 


i 

i 








-23 


45.81 


I 




1 


i 


-11 


45.53 ! 


if 135 ) 
j v o / 

i 


-140 


40.86 j 




+ 8 


45.05 i 


-90 


40.78 




+ 17 


44.93 


1 


-70 


40.67 


i 


+ 27 


44 . 84 


i 


-50 


40.50 


f 


+ 38 


44.72 


i 

i 


-31 


39.60 




+ 47 


44.71 




-20 


38 . 80 


) 

? 


+ 60.5 


44 . 69 


i 


-11 


38.15 




+ 70 


44.69 


! 

i 


- 9 


37.95 


| 




i 


! 


1 


37.15 


( 4 5 V 
V 30 J 




i 


i 


+ 10 


36.90 


-150 


54.67 


l 


+ 27 


36.60 


-130 


54.96 


! 


+ 42 


36.48 




-110 


55.32 


j 


+ 60.5 


36.41 




- 90 


55.46 i 




+ 80 


36.31 




- 70.5 


55.43 ; 


? 








- 50 


55.33 j 


/150) 

V 0 / 








- 30.5 


54 . 84 


-150 


41.29 




- 11 


53.50 


-140 


41.29 

41.21 


+ 8 


52.20 


< 


-100 


+ 27 


51.61 




- 60 


41.08 




+ 46 


51.25 




- 40 


40.79 




+ 60.5 


50 . 88 


i 


- 20 


40.00 




+ 70.5 


50.81 




- 11 


39.54 




+ 81 


50.86 


i 


+ 8 


38.86 




+ 90.5 


50.89 


i 


+ 27 


38.51 








j 


+ 37 


38.48 


( 45, 
V 3 0y 






» 

f i s~ rr \ 


+ 60.5 


38.42 


-150 

-130 


45.35 

45.45 


\(165\ 
j \ 0 ' 


-140 


43.14 


; 

i 


-110 


45.66 


-100 


43.14 


i 


- 90 


45.66 


! 


- 90 


43.12 




- 70.5 


45.62 


j 


- 72 


43.10 




- 50 


45.65 


j 


- 50 


43.00 




- 30.5 


45.34 


I 

i 


- 33 


42.70 




- 11 


45.00 


t 

( 


- 11 


42.15 


\ 

f 


+ 8 


44.40 i 




+ 3 


41.86 > 




+ 27 


44.15 




+ 27 


41.62 ’ 


i 


+ 46 


44.05 




+ 60.5 


41.55 i 


+ 60.5 


43.80 






1 




+ 70.5 


43. 85 



Data as corrected for Temperature recalibration. 
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Field 


Temp . 
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Field 


Temp . 
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(°C) 
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Orien . 


(°cj 
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(* 5 A 
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43.88 


1 (IV 


-150 


41.88 
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43.84 


1 V 4 5 / 

1 


-130 


41.78 


(55. 3\ 

\4 8 . 2 / 




j 


1 


-110 
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-150 
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I 
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41.65 


-140 


55.97 
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41.59 




-130 
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- 50 


41.42 




-120 


56.12 




- 30.5 


41.05 




-110 


56.20 




- 11 


40.35 




-100 


56.34 
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- 90 
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- 80 
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+ 81 
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i 


- 20 
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- 10 
0 
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/ 6 0 \ 
V4 5y 
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+ 20 


54 . 84 
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57.06 | 


1 

l 


+ 40 
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57.00 j 
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- 50 
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1 

1 
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- 30.5 


56.49 


! / 5 5 . 3 \ 
i v-4 8.2y 








- 11 


55.47 


-150 

-140 


42. 37 
42.27 
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+ 18 
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53. 70 


| 


-130 


42 . 22 




+ 27 


53.39 


i 

i 


-120 
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1 

! 
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i 

i 


-100 
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> 

i 
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41.99 




+ 81 
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t 

i 


- 80 
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(90S 

\45J 




j 




- 70 


41.81 


-150 


52 . 32 




- 50 


41.67 


-130 


52.33 


1 


- 40 


41.63 




-110 


52.27 j 


l 


- 30 


41.59 




- 90 


52.22 




- 20 


41.48 




- 70.5 


52.18 




- 10 


41.37 




- 50 


51.99 




0 


40.96 




- 30.5 


51.60 j 
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- 11 
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+ 20 
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+ 8 
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39.90 
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+ 60 
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+ 70.5 


48.47 






j 
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Field Temp. Differences 

Orien. (°C) Taken 



-150 


40.85 


-130 


40 .45 


-110 


40.27 


- 90 


40 . 24 


- 70.5 


40.20 


- 50 


40.14 


- 30.5 


39.39 


- 11 


37.04 


+ 8 


35.51 


+ 27 
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+ 46 
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+ 60.5 
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+ 70.5 


34.55 


+ 81 


34.55 
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34 .54 


-150 


47.14 


-130 


47.01 


-110 


46.86 


- 90 


46.79 


- 70.5 


46.65 


- 50 


46.44 


- 30.5 


46.03 


- 11 


45.38 


+ 8 


44 . 82 


+ 27 


44.47 


+ 46 


44.29 


+ 60.5 


44.16 


+ 70.5 


44.13 


+ 81 


44.05 


+ 90.5 


44.00 
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Coupling Constant (gauss) 




-150 -100 -50 0 50 

Temperature (°C) 



Figure 14a. Experimental results of the coupling 
constant dependence on temperature 
for selected magnetic field orientations, (®) 
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Coupling Constant (gauss) 




Figure 14b. Experimental results of the coupling 

constant dependence on temperature for 
selected magnetic field orientations, (jp 
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Coupling Constant (gauss) 



48 - 




-L 



± 



± 



-150 



-100 -50 0 

Temperature (°C) 
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100 



Figure 14c. Experimental results of the coupling 
constant dependence on temperature 
for selected magnetic field orientations , 
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Coupling Constant (gauss) 




Figure 14d. Experimental results of the coupling 

constant dependence on temperature for 
selected magnetic field orientations, (^). 
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may be due to a crystallographic change in the system occur- 
ring at about 30°C. However, X ray diffraction patterns of 
the crystal compared at approximately liquid nitrogen tempera- 
tures and at temperatures above 30 °C showed that no such phase 
change had occurred. 

The experimental data shows that for those magnetic field 
orientations at which the two magnetically inequivalent sites 
can be distinguished there is a characteristic rise or drop 
in the coupling constants after leveling off as the tempera- 
ture is lowered. This drift of the coupling constants remains 
unexplained. The total shift in coupling constants tabulated 
in Table II are based on the shift observed between the points 
where leveling occurs at low and high temperatures , the noted 
drift being disregarded. This procedure proved to be -justi- 
fied based on calculations of the overall shift to be discus- 
sed under the relaxation effect results. 

B . MODEL DEVELOPMENT 

Experimental proof of the coupling constant dependence 
upon temperature led to development of a model by which the 
effect could be explained. The work of Tolies et.al [Ref. 

21] and Ohigashi and Kurita [Ref. 15] showed rotation of the 
CH 2 group about the C-C bond; hence the restricted rotor was 
chosen as the initial model to be studied. 

The Hamiltonian for the restricted rotor has been discus- 
sed in Section III and was seen to be 

H = BL^ + j(l-cos2<]> ' ) 
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TABLE II. Total coupling constant shifts observed experi- 



mentally 



Field Orientation 

! 


] 1 

Total Shift 
(gauss) 


j 

( 9 CL 
V 0 )' 

j 


! ; 

0.35 

t : 

i i 


(120\ 

Vo/ 


1.25 ; 

i I 


f 1 * 5 ) 1 

V o * 


4.48 : 

« 


{ 150' 
0 ■ 


2.90 


(165': 

v 0 / 


[ 

1 1.60 


O 03 


5 

1 .86 
i 

1 


i M5\ 

V3oy 


. 

4.61 


/' 45 > 
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! 1.80 ! 

i ; 


/ 5 5 . 3\ 
V4 8 . 2) 


J i 

3.95 ; 

j ! 


{ 55. 3N 
! Y4 8 . 2/ 

j 


♦ 

2.34 


(S?) 


4.28 j 


<£l) 


( i 

2.50 ; 

1 1 


(90\ 

V4 5/ 


3. 70 j 

i i 

1 > 
i 




| 5.69 

! ! 


(» 


1 

2.70 1 

\ ) 

i 



^"Measured by the difference in extreme absorption lines. 
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where solution of the Hamiltonian matrix elements and comout- 



erized di agonalization of the resulting secular determinant 
led to the eigenvalues characteristic of the restricted rotor 
system. Using the results of Tolies et.al. [Ref. 21] for an 
H-C-H angle of 116.9° and a C-H bond length of 1.08A, the 
moment of inertia, I, of the CH 2 rotation about the C-C bond 
was calculated to be 1.69408A. The barrier ootential, V, 
of 5.8 kcal/mol was converted to units of V/H and found to 
be 60.8524 THz (terahertz). Knowing I, B was calculated to 
be 0.298414 THz in units of B/h . From these values, the 
Hamiltonian matrix elements were readily calculated for 
integer values of m. 

As it was obviously impossible to work with a determinant 
of infinite size corresponding to the range of m values of 
minus to plus infinity, several determinants of increasingly 
larger size were solved until the variance in eigenvalues 
from one size to the next larger size was within 0.2% or 
approximately 0.1 THz. It was found that a 29 x 29 deter- 
minant (m having values -14 to +14) was the smallest to meet 
these requirements. Tabulation of the eigenvalues for the 29 
x 29 determinant appear in Table III. 

Examination of the resulting eigenvalues for the restrict- 
ed rotor shown in Table III showed each energv level to be 
doubly degenerate as expected from the solution of Schroed- 
inger's equation which yields energy values 

_ m 2 h 2 
E - ~~2I 

so that a -m or +m value yields the same eigenvalue. 
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TABLE III 



Eigenvalues of the restricted rotor model . 



Energy | 
Level 

i 

i 


| Restricted Rotor 

Eigenvalues 1 
THz) 


! 

0 


4.1853 


[ i 


4.1853 

i 


i 

i 


12.3979 | 


\ 

\ 

1 


12 . 3981 

i 


2 


20.2843 | 


j 


20.2861 


i 1 

i 3 

1 i 


27. 8221 

27.8318 i 

. i 


i ; 

1 4 

; 1 

! 

j 


34.9854 

35.0223 

! 


I 5 

t J i 

i 


41.7450 

41.8451 


i | 

| 6 
j 


48.0342 

48.2856 j 


7 


53.9234 

54.0146 

t 


8 


58.1838 ! 

60.0289 j 

1 



1 



Values obtained using a barrier to internal rotation 
5 . 8 kcal/mol . 



of 
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Further examination showed that for the first four levels 
above ground state there was an approximately 8 THz difference 
between levels and that the ground state level itself was 
4.19 THz which is approximately half of the energy level dif- 
ferences . This observation suggested immediately that the 
restricted rotor system might well be approximated by an 
harmonic oscillator which has energy levels represented by 

E = (n+i-)hv (equation III-A2-1) 

and which in itself is one limit of a restricted rotor. The 
harmonic oscillator approximation was applied to the restrict- 
ed rotor, and energy levels were calculated by using equation 
III-A2-1 as rearranged to obtain the eigenvalues in units 
consistent with those of the restricted rotor. (Equation 
III-A2-3, Page 27 ) . Values of n equal zero to eight were 
used in the calculation, the results of which were compared 
to the restricted rotor results and appear in Table IV. Ex- 
amination of Table IV shows that there is good agreement be- 
tween the two models in the first four levels above ground 
state, but that at higher energy levels, the oscillator model 
yields increasingly higher results. 

A Boltzman distribution calculation to obtain approximate 
population of energy levels in the harmonic oscillator is 
shown in Table V. 

Based on negligible population of energy levels higher 
than the fourth excited state, the harmonic oscillator was 
assumed to be a good approximation to the restricted rotor 
and subsequent calculations were based on this approximation. 
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TABLE IV. Eigenvalues of the restricted rotor and the 



harmonic oscillator model 



Energy 

Level 


Restricted Rotor 
Eigenvalues 
(THz) 


Harmonic Oscillator 
Eigenvalues 

(THz) i 


0 


4.1853 

4.1853 


! 

| 

4.2614 

< 


1 


12. 3979 
12 . 3981 


12.7842 


2 


20.2843 
; 20.2861 


21.3070 j 


3 


27. 8221 
27.8318 


29 . 8298 

{ 


t 

[ 4 
1 


34.9854 

35.0223 


38.3526 

1 


[ 5 

i 

? 

i 


41.7450 

41.8451 


46.8574 ; 

i 

1 


3 

6 


48.0342 

48.2856 


55.3982 

j 


7 


53.9234 

54.0146 


63.9210 

1 

j 

i 


8 


58.1838 

60.0289 


72.4438 j 

1 

i 



1. Values obtained using a barrier to internal 
rotation of 5.8 kcal/mol. 
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TABLE V. Energy level populations in the harmonic oscillator 



model . 



Energy Level 

| 


Percentage 
Population 
At 273°K 


i j 

Percentage j 
Population I 
At 133°K J 


0 


68.9 


94.54 


1 ; 


1 i 

23.5 

! ! 


5.15 

[ 


2 ! 


1 5.8 


I 0.29 

| j 


3 


1.6 


i 0.02 1 

i 


4 

) 


0.2 


0.00043 

I 

! 



C. SHIFT OF COUPLING CONSTANTS DUE TO DIPOLAR INTERACTION 



In Section III development of the calculation used to 
determine the average coupling constant as a function of 
temperature (<A> T ) based on the effect of dipolar interaction 
was discussed. It was seen that the relation to be solved 
was 



A- +9 0 



<a> t - \ 



-90 



3.77629A , £ e H 



-E /KT 
n' 



n 



nU)' 



-d<|)' 



(equation III A- 3-27) 

In the development of this equation, a number of comouterized 
subroutines were written to the basic computer orogram SSESR 
to yield program PREDICT II which appears in Appendix A. Not 
considering hyperconjugation effects (Subroutine CCCA) which 
will be considered in the next section of this work, equation 
III-A3-27 was solved by PREDICT II to vield the shift in 
coupling constants for selected magnetic field orientations 
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over a temperature range of +90°C to -150°C. The overall 
shifts are tabulated in Table VI. 

Before considering the significance of the results, the 
manner in which PREDICT II was utilized will be discussed. 

The theoretical development of each of the individual 
subroutines was discussed in Section III-A-3. Basic input 
data to the main program consisted of those parameters which 
Tolies et.al [Ref. 21] determined for the planar form of the 
Ct^CC^ radical plus parameters determined in this work to 
include , 

The number of protons in the radical = 2. 

The magnitude of the magnetic field or 
the nuclear Zeeman frequency = 13.8 THz. 

The principal values of the coupling tensors (coupling 
constants) as listed on page 18 of this work. 

The orientation of external magnetic field, (®) . 

The direction cosines of the H tensor as listed on 
page 20 of this work. 

The first five eigenvalues of the harmonic oscillator 
as developed in section V-2 and appearing in Table IV. 
The temperatures at which the coupling constants were 
to be calculated. 

The limits of integration, -90° to +90°. 

The increments of d<J) ' used in the numerical integra- 
tion, 3° . 

Miscellaneous constants as appropriate. 
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The output yielded the esr spectrum at any given angle of 
rotation, <j> ' , in terms of transition frequencies and inten- 
sities as intermediate results. The coupling constants were 
then determined by frequency differences and automatically 
fed back into the program to yield the average counling con- 
stants <A>,p at given values of input temperatures for dioolar 
coupling effects. 

It was of prime concern to insure that PREDICT II was 
rotating (operating) properly and yielding reliable results. 
Hence, two operational tests were devised and used as checks. 

1 . Test One for Proper Rotation 

The spectra computed at +90° rotation miqht be ex- 
pected to be coincident with that at -90° , and this would 
make an ideal check for proper rotation. Initially the re- 
sults of the two rotations did not agree which led to 
incorporation of subroutine REFH which insured orthonormali tv 
of the coupling tensors and H 2 . It must be noted, in 
addition, however, that the input coupling constants for pro- 
ton one are not exactly the same as those for proton two, 
hence the dipolar effect cannot be expected to yield iden- 
tical results for +90° and -90° of rotation. To check for 
proper rotation, both sets of input coupling constants were 
made identical (arbitrarily, identical to those of proton one) 
and rotation by PREDICT II was found to yield identical spec- 
tra at +90° and -90° of rotation. This procedure gave one 
good confirmation that PREDICT II was operating properly. 

For all future calculations except where indicated, input 
coupling constants were assigned their original values. 
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2 . 



Test Two for Proper Rotation 



Using the transformation 

; cos<J> sin<}) 0 \ 

( i ' . - 1 

! M ; -simj) cost}) 0 ; :M 

\ ' j i 

i o o i / 

180 

an orientation of magnetic field ( ^ ) was rotated through 

an angle cp of 90° , where the matrix M is that for the direc- 

18 0 

tion cosines of magnetic field ( q )• This transformation 
was performed with aid of a desk calculator and yielded new 
direction cosines for the magnetic field of 

0.323832 
0 . 896897 
0.301199 . 

These direction cosines were then fed to PREDICT II as input 
direction cosines for the magnetic field while also setting 
the input coupling constants identical as in test one, and 
the spectra were computed for angles of CH„ rotation of -90° 
and +90°. If the program was rotating properly, the result- 
ing spectra were expected to correspond to that spectrum for 

180 

an input of magnetic field orientation ( q ) at 0° rotation 
of the CH 2 group. Agreement was excellent and it was conclud- 
ed that PREDICT II was rotating properly. 

Table VI shows that a shift in coupling constants be- 
tween -150°C and +90°C did occur due to dipolar coupling, 
however the magnitude of the effect in no instance was that 
observed experimentally. The direction of the dipolar cou- 
pling effect was seen to be dependent on the external magnetic 
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TABLE VI. Calculated coupling constant shift over tempera 
ture range of -150°C to +90°C due to dipolar 
coupling. 1 



Field Orientation 


Shift (gauss) 


(30) 


; 

-0.013 ! 

i 


© 


s 

+0.006 i 

' l 

1 ► 


CTi G\ 
O O 


1 

0 .000 


f 75X | ' 


+0 .033 

i 

\ 


v 3 0 / 


( 75 i 

\6oy 

j \ 


-0.035 j 

\ 


: 90 ) 

UsJ ; 


j 

-0.054 

j 


/ 9 ON i 

V-4 5y 


+0.045 

! 


( 9 °) 
! \75y 


! I 

-0.091 S 

i 


t ( 90\ | 

V9oy 


! 

-0.079 

l 

1 


\ 

e?°) 


+0.048 


Co 3 !) 


+0.013 


(T) 


1 

-0.013 

! 


er) 


i 

-0.079 

i 


Co 0 ) 


-0.119 



"''All shifts are determined by subtracting the coupling con 
stant at -150°C from that at +90°C. 
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field. Table VI shows that some magnetic fields caused an 
effect which led to an increase of coupling constants with 
increasing temperature (dipolar coupling values positive) 
while others led to a decrease in coupling constant with in- 
creasing temperature (dipolar coupling values negative) , the 
general effect observed experimentally. This effect was also 
observed graphically with plots as shown in Figure 15 where 
the calculated coupling constants due to dipolar interaction 
have been plotted against angle of CH ^ rotation, <j> ' , for 
several magnetic field orientations. 

As the temperature is raised, greater angles of b ' 
are attained. Hence, analysis of Figure 15 showed that in- 
creasing temperature would be expected to cause an increase 

in coupling constants for those plots having positive curva- 

120 

ture at <j) ' equal zero (e.g. the ( ^ ) plot) and a decrease 

in coupling constants for those plots with negative curvature 

90 

at b ' equal zero (e.g. the ( 75 ) plot) . The curvature of the 

plots and the direction of coupling constant shift in Table 

150 

VI agreed. The ( Q ) orientation showed curvature of nearly 
zero and this too agreed with the nearly zero value seen in 
Table VI. 

These results showed that although some shift in cou- 
pling constants occurred due to dipolar interaction, the 
direction of shift was not always consistent with that ob- 
served experimentally and hence dipolar interaction was not 
the primary cause of the effect. 
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Figure 15. Calculated coupling constants due to dipolar interaction as 



It was of note that although the calculated shift was 
small, the magnitude was the same order as that observed for 
the methyl radical by Zlochower et.al [Ref. 23]. 

D. SHIFT OF COUPLING CONSTANTS DUE TO FERMI CONTACT INTER- 
ACTION 

The theory of hyperconjugation and resulting Fermi contact 
interaction was discussed in Section III-B. 

Using a computerized orbital calculation based on the CNDO 
approximation [Ref. 19] the proton 2S molecular orbital coef- 
ficients were calculated for various angles of rotation, 4 > ' . 
The coefficients squared yielded the amount of free electron 
spin density that was imparted to the proton 2S orbitals for 
any angle of rotation. The solid curve of Figure 16 shows a 
plot of the resulting spin densities versus the angle of rota- 
tion, <f> ' , which upon examination appeared to closely approx- 
imate a function represented by 

= Asin 2 2cf> ' 

-2 

of amplitude A = 0.273 x 10 . The dotted curve of Figure 16 

is a plot of this function to show that, in fact, the hvoer- 
conjugation effect did very closely approximate the sin 2 2<j)' 
relationship. 

If the entire electron (spin density = 1.0) were found in 
a proton 2S orbital, then the coupling constants of the planar 
radical (where due to orbital symmetry, no spin density ap- 
pears in the proton 2s orbital) would be increased bv the 
factor B = 1,420.4058 MHz, the hyperfine coupling constant of 
a hydrogen atom with a complete electron in the 2s orbital. 
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Figure 16. Results of molecular orbital calculation for Fermi contact 
(hyperconjugation) interaction compared to f , , . Asin 2 2(J) ' . 



Hence, the input coupling constants for the planar radical, 

, were adjusted for Fermi contact interaction throuoh the 
relation 

A., = A. + B(0.273xl0 _2 sin 2 2<j) 1 ) . 

4> i 

This relation was incorporated into PREDICT II through sub- 
routine CCCA (Appendix A) and values of the average coupling 
constants <A>^, were again calculated having both dipolar and 
Fermi contact interaction incorporated. Column 2 of Table VII 
shows the results obtained for the total shift in coupling 
constants over temperature range -150°C to +90°C which can be 
compared with column 3, the values for dipolar interaction 
only . 

Since the hyperconjugation effect has no dependence on 
orientation of magnetic field but only on angle of rotation, 

4> 1 , it was expected that the effect should be constant for 
all field orientations. By subtracting the dipolar effect 
from the total effect, the magnitude of hyperconjugation 
effect on the coupling constant shift was determined, the re- 
sults being tabulated in column 4 of Table VII. As expected, 
the effect is essentially constant for all magnetic field 
orientations having a value equal to about -0.19 gauss over a 
temperature range of -150°C to +90°C. The negative value 
corresponds to a decrease in coupling constant with increas- 
ing temperature. 

The calculated average coupling constants incorporating 
both effects were plotted versus temperature giving plots 
characteristic of Figure 17 for all magnetic field orientations 



90 



TABLE VII. 



Calculated hyperconjugation effect (Fermi con- 
tact interaction) on coupling constant shift. 1 



Field 

(orientation 


Shift with 
Dipolar and 
Fermi Contact 
Interaction 
(gauss) 


Shift with 
Dipolar 
Interaction 
Only 
(gauss ) 


Shift Due to 
Hyperconj ugation 
(Fermi Contact) 

(gauss) 


(8) 


-0.209 


-0.013 


-0.196 


Go) ! 


-0.189 


-0.006 : 

! 


-0.195 

I 


© 

l 


-0.196 


0.00 


I 

-0.196 


^ 7 5 \ 

Voo; 


-0.158 


+0.033 


-0.191 


; Go) j 


-0.220 


i 

-0.035 


. 

| -0.185 


! («) j 


i 

-0.245 


-0.054 

< 


-0.191 

i 


( 90\ 

i V-4 5 ) 

- > 


! i 

-0.128 


i 

+0.045 i 

i 


-0.173 


( 9 °\ 

| V75/ 


-0.283 


-0 . 091 


; | 

-0.192 

i j 


/90N 

\90j 


i 

-0.264 


-0.079 


! -0.185 

1 

( 


(G°) 


-0.128 


1 

+0.048 


t j 

-0.176 

j 


O U> 


-0.167 


+0.013 

| 


: 1 

-0.180 i 

1 


(G°) 


-0.198 


-0.013 j 


i 

-0.185 j 


CG 5 ) 


-0.268 


1 

-0.079 


1 1 

1 

-0.189 


O 00 


-0.311 


+0.119 

! 


j 

-0.192 

! 1 



'*'A11 shifts are determined by subtracting the coupling 
constant at -150°C from that at +90°C. 
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studied. Examination of these plots and the data of Table 
VII led to several conclusions regarding the model as calcu- 
lated to this point. 

First there appeared to be a definite effect on the cou- 
pling constants when both dipolar and Fermi contact interac- 
tion were considered. The overall effect led to an increase 
in coupling constants with decrease in temperature as observed 
experimentally, the primary effect arising from Fermi contact 
interaction. The plots of Figure 17 indicate a decrease in 
the temperature coefficient as lower temperatures were ao- 
proached indicative of quenching of the rotational motion. 

This observation agreed with experiment but it was of note 
that no such temperature coefficient decrease occurred at 
higher temperatures in the calculated effects. Of greatest 
significance, however, was the magnitude of the calculated 
effects which was seen to be only about 0.1 to 0.3 qauss, 
not nearly that observed experimentally (1-6 gauss). 

Figure 18 compares a typical experimental result with that 
calculated on the basis of dipolar and Fermi contact 
interaction . 

Based upon the foregoing observations , it was concluded 
that although the calculated model showed an overall effect 
generally of the form observed experimentally, and althouqh 
dipolar and Fermi contact interaction can be expected to con- 
tribute to the shift in coupling constants, these were onlv 
a small part of some larger effect which was yet to be 
established. 
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Coupling Constant (gauss) 




Figure 17. Dependence of coupling constants on tempera- 
ture due to dipolar and hyperconjugation 
effects calculated. 
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Coupling Constant (gauss) 




Temperature (°C) 



Figure 18. Calculated dependence of coupling constants 
on temperature due to dipolar and hyper- 
conjugation effects compared to experimen- 
tally observed temperature dependence. 
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E. SHIFT OF COUPLING CONSTANTS DUE TO FLAP MOTION 

Out-of-plane flapping motion was discussed in Section 
III-C and was seen to be a type of hypercon j uqation effect 
leading to Fermi contact interaction as with the rotational 
model but under a different mathematical relationship. An 
attempt was made to determine if this type vibration was 
occurring in the CH 2 CO 2 radical and the magnitude thereof. 

A molecular orbital calculation of the same type used 
for the rotational model [Ref. 19] but with various angles 
of flap led to the proton 2S orbital spin densities tabula- 
ted in Table VIII. In addition, the relative energies of 
the radical at the various angles of flap were also deter- 
mined from the molecular orbital calculation and appear in 
Table VIII. 

No simple function relating the coupling constants to 
angle of flap could be obtained from the data of Table VIII. 
However, a plot of the energy versus angle of flap as seen 
in Figure 19 showed a double minimum which indicated that in- 
deed, a flapping vibrational type of motion may exist within 
the system. The fact that the energy of the radical de- 
creases to a minimum at about ±20° of flap indicated that if 
flap does exist, then the most stable form of the radical may 
not be planar as reported by other workers. [Refs. 15 and 21]. 

A means by which this possibility could be tested would 
be to lower the temperature of the crystal to liquid helium 
values in an attempt to freeze out the two stable forms , one 
with an angle of flap of about + 20° the other of - 20° . An 
esr spectrum would then be expected to disclose four 
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TABLE VIII. Results of molecular orbital calculations for 
flap vibrational motion. 



Angle of Flap"*" 
(degrees) 


2 

H„ Spin Density 
/s x lO- 3 


Relative Energy 
of Radical (aeu) 


0 


0.000 


-52.10978 


1 


0.004 


-52.10980 


3 


0.049 


-52.10990 


5 


0.090 


-52.11006 


10 


0.182 


-52.11058 


15 


0.181 


-52.11097 


20 


0.133 


-52.11110 


22 


0.118 


-52.11106 


24 


0 . 096 


-52.11097 


26 


0.071 


-52.11081 


28 


0.054 


-52.11059 | 

t 


30 


0.040 


1 

-52.11030 


32 


0.030 


1 

-52.10994 


40 


0.003 


-52.10773 ! 

1 


50 


0.001 


-52.10301 


60 


0.002 


-52.09578 


70 


0.000 


-52.08570 


80 


0 . 020 


-52.07244 


90 


0.129 


-52.05565 



Negative angles of flap yield nearly identical results. 
Values reported to three significant figures. 
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Figure 19. Relative energy of the Cf^CC^ radical at 
configurations caused by CI^ group out- 
of-plane flapping motion. 
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inequivalent sites rather than the two observed at the temper- 
atures studied in this work. If the four sites were to be 
identified, then raising the temperature through increments 
would yield data from which, via Arrhenius plots, a barrier 
to internal flapping motion could be determined. 

Flap vibration is not the most likelv mechanism to ex- 
pect. No experimental evidence of the effect was observed 
over the temperature range of -150°C to +100°C. Based on the 
more obvious rotational mechanism which can be substantiated 
by experimental data, it was concluded that if flapping does 
occur, it is most probably characterized by a very low fre- 
quency and large amplitude with insignificant effect upon 
the system's coupling constants. 

F. SHIFT OF COUPLING CONSTANTS DUE TO RELAXATION EFFECTS 

The development of the relaxation model and methods of 
calculation were discussed in Section III-D. It should be 
noted that the development was based again on the original 
model of the C^ group rotating about the C-C bond of the 
CH 2 CO 2 radical. 

The first step in making the calculation was determina- 
tion of the limiting conditions and hence the overall shift 
in coupling constants between these conditions. SSESR with 
a slight modification was used to calculate the limits. 

For a given magnetic field orientation, the relative 
coupling constants at the low temperature limit (low flip 
probability) were obtained through solution of Box's [Ref. 2] 
relation. 
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the results of which appear in Table IX column 2. 

The high temperature limit was based on the approximation 
that at high temperature conditions , the rate of flip is so 
rapid that the effective field experienced by each proton will 
be an average of the two fields h^ and as expressed bv the 



relation 



_|h l( + ) +h 2( + ) 1 + 



h l (-) +h 2 (-) 



Results from solution of this relation are seen in column 3 
of Table IX. 

The difference between the two limiting coupling constants 
yielded the overall shift A c (calculated shift). Since the 
program calculated the values in units of megahertz it was 
necessary to divide the resulting overall shift by the con- 
version factor 

2 . 80 MHz/gauss 

in order to compare the results with experimental overall 

shifts (A ) . Column 4 of Table IX shows the final calculated 
e 

shift and column 5 the overall shift observed experimentallv . 

The values of A g were plotted versus the values of A c> 

As seen in Figure 20 , the points cluster about a straight 
line, deviations being assessed to other effects which have 
not been isolated. 

Each point on Figure 20 represents a different magnetic 
field as tabulated in Table IX. It is of note that approx- 
imately ten hours of experimental work and computations were 
required to establish one point. It was noted, in Section V-A 
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TABLE IX. Comparison of overall calculated coupling constant 
shift with experimentally observed shift. 



Field 


Low T Limit 1 


High T Limit 1 


A 

c 


A 

e 


Orien. 


(MHz) 


(MHz) 


(gauss) 


(gauss ) 


to 2 


113.4231 


111.3107 


0.745 


0.40 


to) 


113.0746 


98.9371 


5.049 


4.25 


TO) 


113.9912 


101.2032 


4.567 


4.48 


TO) 


116.3903 


180.6752 


2.755 


2.90 


TO) 


122.9638 


118.3930 


1.632 


1.60 


TO) 3 


132.6437 


127.5343 


1.825 


1.86 


($ 


157.8302 


146.4922 


4.049 


4.61 


TO) 

A 


127.5728 


122.4376 


1.834 


1.80 


( ' 55 - 3 v 

V 48 . 2 y 


161.8221 


149.7523 


4.311 


3.95 


4 


118.4243 


112.6443 


2.064 


2.34 


GO 


161.2777 


149.2828 


4.284 


4.28 


TO) 


117.4068 


110.2377 


2.560 


2.50 


01) 


147.5643 


138.9060 


3.092 


3.70 


TO) 


113.6078 


96.3793 


6.153 


5.69 


oil 


133.2189 


128.0851 


1.834 


2.70 



Calculated limiting conditions 
Orientation of magnetic field along 
Orientation of magnetic field along 
Orientation of magnetic field along 
Orientation of magnetic field along 



the c crystal axis 
the a* crystal axis 
the C-C bond 
the b crystal axis 
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that experimental plots of coupling constant versus tempera- 
ture for those cases where <J> yL 0 (®) , showed drift in the 

coupling constants with lower temperatures after they leveled 

off at about -60°C. The values of A used in Table IX and 

e 

Figure 20 were, for these cases, measured between the high 
and low leveling off points, and as seen in Figure 20 agree 
well with the calculated shifts. Initially, the deviations 
from linearity of Figure 20 were thought perhaps to be due to 
the drift observed experimentally. However, examination of 
the data showed no general correlation to confirm this fact. 

Based on the linearity of Figure 20 it was felt that the 
method of calculation for the limiting conditions was repre- 
sentative of the system. Having determined the limiting 
conditions and hence the overall shift in coupling constants, 
the calculation of coupling constants for intermediate flip 
probability (intermediate temperature) was undertaken using 
the Monte Carlo technique discussed in Section III-D. 

The Monte Carlo technique is based upon random statistics , 
hence the more points on a given curve that are calculated 
and the more times each point is calculated, the more reliable 
the final results are expected to be. Numerous calculations 
of the damped cosine function (Equation III-D4-1, Page 52 ) 
approximating the M component of magnetic moment and the 
resulting frequency were made by varying the input parameters 
of flip probability, d, (from 2 to 100), precession increment, 
A, (from 2 to 10) , number of cycles calculated (1 or 2) , ran- 
dom number input, and the number of times a particular set of 
parameters was recalculated (100 to 400 times). The angle a 
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A c (calculated) (gauss) 



Figure 20. Comparison of total coupling constant shift 

calculated (A ) with that observed exoeri- 
c 

mentally (A g ) . 
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between the Z axis and the effective magnetic field was chosen 
as 25.84° which permitted frequency to be represented as rela- 
tive values between 0.90 and 1.00 for a given flip probabil- 
ity. The magnitude of the results were independent of a, but 
the specific value was chosen for ease in converting the fre- 
quency data to increments of overall shift in the coupling 
constants . 

For each value of d and A used to determine the best fit- 
ting damped cosine function, values of correlation time, t , 
were calculated by equation III-D6-4 based on a frequency of 
precession, v, of one hertz. The values of t c calculated 
thusly were plotted against the frequency of the best fitting 
damped cosine curve with the results shown in Figure 21. 

There is a considerable scattering of points in Figure 21 
due to the random statistical method of calculation. However, 
it was not difficult to describe a curve most representative 
of the greatest majority of points. The "S" shaped curve 
was noted to be most characteristic of that type function ob- 
served experimentally. 

For further calculations , only those points lying on or 
close to the curve described in Figure 21 were used, these 
conditions being considered as most representative of the re- 
laxation model. Examination of the spectra calculated by 
PREDICT II for given temperatures and magnetic field orienta- 
tion showed that the average frequency of precession eXDeri- 
enced by the system was about 30.28 MHz. Hence, each 
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Figure 21. Relative frequency as a function of correlation time based on 
correlation time calculated on a 1 Hz basis. 
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correlation time as calculated based on one hertz was adjusted 



to a 



real 



consistent with 



r i t 
real 



the rea 

/30. 28 
c 



1 

X 



system where 
10 6 . 



Using equation III-D6-5 (Page 57 ) and ' T rea ]_/ values of 
temperature were calculated to correspond to the frequency 
for a given probability d, the results of which are tabulated 
in Table X. 

The plot of relative frequency versus temperature in Fig- 
ure 22 shows that the "S" shaped function characteristic of 
the observed experimental function was retained. 

The final step was to convert the data of Figure 22 or 
Table X to coupling constant shift for specific maqnetic 
field orientations. As seen in Figure 22 the frequency was 
represented as a relative frequency shift between 0.90 and 
1.00 where at any given temperature the percentage of total 
shift between the limiting frequencies and hence couoling 
constants could be readily determined. For example, at a 
temperature of 0°C, 43% (0.943) of the total shift in cou- 

pling constants is expected to have taken place. Since verv 
good agreement was obtained between the experimental and 
calculated low temperature limits and overall shift of cou- 
pling constants, the percentage shift at a given temperature 
was applied to observed overall shifts and curves were plot- 
ted based on the calculations just discussed, the results 
being compared with experimental curves as seen in Fiqures 
23a through 23c. 
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TABLE X. Results of Monte Carlo calculations correlated 
to temperature. 



[Probability 
Factor (d) 


Relative Frequency 
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Temperature (°C) 
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Relative Frequency (Hz) 



1.00 




Temperature (°C) 



Figure 22. Relative frequency as a function of 
temperature based on an average 
precession frequency of 30.28 MHz. 
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Figure 23a. Comparison of calculated dependence of 
coupling constants upon temperature 
due to relaxation effects with 
experimentally observed dependence for 

16 5 

a magnetic field of ( q ) . 

I Experimental; O Calculated. 
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Figure 23c. Magnetic Field (qq) . 
I Experimental; O Calculated. 
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Comparison of the calculated and experimental curves was 
made by noting the temperatures at which the temperature coef- 
ficients changed sharply (temperatures where the curves level 
off) and the inflection points. 

The calculated curves followed a function closely approx- 
imating the experimental curves, but with a shift upward in 
temperature of 10°C to 20°C. The position of the calculated 
curves depended heavily on the value of T rea ^/ shifting down- 
ward in temperature for lower values of average frequency 
experienced by the system. However, this fact was of minimal 
value in that from calculated spectra results, the average 
frequency of 30.28 MHz was the proper value to use for the 
system. 

Although the calculated curves did not agree with experi- 
mental curves as well as would have been desired, based on 
the random statistical method used in making the calculation, 
it was felt that the results obtained were in fact very good 
in that it was shown that the relation observed could be 
essentially constructed entirely from theoretical calculations. 

Of greatest significance, however, was the fact that the 
temperature dependence observed in the CH 2 CO 2 radical esr 
spectrum could now be explained as being due primarily to a 
relaxation effect as was developed in this work. 
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VI. CONCLUSIONS 



As originally noted by Crawford [Ref. 4] a marked temper- 
ature effect upon the esr spectrum of the CI^CC^ radical 
formed in X-irradiated zinc acetate dihydrate was verified. 
Decrease in temperature corresponded to increase of coupling 
constants extracted directly from the spectra of a given mag- 
netic field orientation. The total shift in coupling con- 
stants does not spread over the entire temperature range of 
-150° to +90° studied but was observed to occur between the 
rather limited range of -60°C to +30°C. 

Development of a mathematical model to describe the ob- 
served effect led to considerations of several phenomena 
expected to take place in molecular systems, the results of 
which indicated that the effect was due to a combination of 
factors . 

Dipolar interaction due to rotation of the group about 

the C-C bond agrees closely with that observed by other work- 
ers for a rotating methyl group but does not account for the 
observed experimental effect. Fermi contact interaction due 
to hyperconjugation effects in the rotational model shows cou- 
pling constant shifts with temperature of the general nature 
observed experimentally, but again, the magnitude of the effect 
accounts for only a very small portion of the observed shift. 

Flap vibrational motion, although perhaps taking place in 
the system as indicated by molecular orbital calculations, is 
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not observed experimentally and is not expected to be a fac- 
tor of any great importance on the coupling constants. 

The primary effect causing the observed temperature depen- 
dence is due to a relaxation mechanism within the molecular 
system. This effect is confirmed experimentally and a model 
has been developed using random statistical calculations. 

It is felt at this time, that given an orientation of exter- 
nal magnetic field, and a temperature at which the CH^CO^ 
radical is stable, the coupling constant can be calculated 
to within ±1.0 gauss. 

The procedures used to studv the temperature effect are 
not restricted solely to the zinc acetate system but could be 
applied to other systems in which similar relaxation proc- 
esses are expected to occur. 

Although it appears that this work combined with the work 
of Tolies et.al [Ref. 21] and Ohigashi and Kurita [Ref. 15] 
nearly exhausts the Ct^CC^ radical system, there are yet a 
few remaining studies which would be worthwhile to undertake. 
Flapping motion has not been confirmed nor proved to be non- 
existant. The drift in coupling constants at low temperatures 
in some orientations of magnetic field has not been explained. 
Each of these two factors deserves an explanation but of most 
value to this work as well as to a general contribution to 
the field of molecular studies would be the development of a 
rigorous mathematical model by which the coupling constants 
could be calculated over the entire range of temperature. 
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This is a challenge yet to be met in the field and it is the 
hope of this author, that this work may be a major step in 
that direction. 
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APPENDIX A 



A computer program, PREDICT II, for calculating the av- 
erage coupling constants of the radical as a function 

of temperature. Various combinations of the subroutines may 
be used to calculate intermediate results as was discussed 
in the foregoing thesis. 
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SUBROUTINE AVBZMIAAVGTlt AAVGT2, EN,T»FAVe) 
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RETURN 

FNO 



A sample of output for PREDICT II 



NUCLEAR ZEEMAN FREQUENCY = 13.80000 



COUPLING CONSTANTS 

-26.80000 -59.59000 -92.15000 

-26.41000 -58.98000 -92.78000 



THETA AND PHI 

45.0000000 30.0000000 



DIRECTION COSINES FOR H 

-0.7071068 -0.3535534 0.6123724 

TEMPERATURE 

363.0000000 



AVG COUP CONST (1) 


, EN 


LEVELS 1-5 








0. 261480556D 


02 


0 .259451835D 


02 


0.2581 36 737D 


02 


0.257394 9 31D 


02 


0.257109 819D 


02 






AV COUP CONST (2) , 


EN 


LEVELS 1-5 








0.29 804999 3D 


02 


0.294510057D 


02 


0 . 291776688D 


02 


0. 289694564D 


02 


0 .288132342D 


02 







average coupling constants 

0 . 260673688D 02 



0.296566194D 02 
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APPENDIX B 



A computer program to calculate spin relaxation effects 
using a Monte Carlo technique. 
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CCMPUTFf. PROGRAM RpR V CNTE CARLO CALCULAT I CN 



C NP IS MR PF POINTS CALCULATED 
C NT IS NP °F T I w 'f S EACH POINT IS CALCULATED 
C IX IS PAN-DON N U M P E 0 INPUT 
r IP IS PRO B API! I T Y 

C AL DPSCRIPES ANGLE ALPHA CF EFFECTIVE FIELD 
C DEL IS I NO RE N’T NT CF PRECESSION 

C TH AND PH APE ANGLES GIVING INITIAL MOMENT PCSITICN 
C 



1 

100 
10 CO 

G 



102 
1 C 3 



20 

21 



ECC> 

RM ( 3 , 3 > 



D I HE MS I PM A ( 500 , IA( 

DIMENSION X( 3) , XP( 3 ) , 

READ 100, NP,NT,IX,IP,AL,DEL,TH,PH 
FORMAT ( 215,21 1C,4F5.2> 

IF(NF) 1CCC»10CC,5 
STOP 

CONTINUE 

W C ITE(8*1C3) NP,NT,IX,IP,AL,36L,TH,PH 

TP = IRC. 0/3. 1415R27 

AL = A I /TP 

DEL. = PEI /TP 

Th « TH/TP 

PH = PH/TR 

TH S = TH 

PKS = PH 

CC = CnS(DEL) 

SD = SIN(OEL) 

CA = COS(AL) 

SA = SI N( At ) 

PM 1 , 1 ) = CD 
R N ( 1 » 2 ) = 

RH( 1 ,3 ) = 

PH(2,n = 

R N ( 2 , 2 ) = 

PH (2,3) = 

R N ( 3 , 1 ) = 

PH (3,2) = 

R N ( 2 * 3 ) = 

rr io i = 

A ( 1 ) = O.C 

PRINT 1C2 , NP , NT, I X , I F , AL ,DEL ,TH, PH 
FORMAT ( IX, 215, 21 10,4F5.2) 

FORMAT ( // , IX ,215 ,2 I 10 ,4F5. 2 ) 

DP 51 IC = 1 , NT 
X ( 1 ) = SI N( THS > *00$ ( PHS ) 

X < 2 ) = SI N< THS )*SIN< PHS ) 

X( 3) = CCS( Tw$ ) 

rr 50 ja = i , n p 
nr 20 i = l ,3 

X F ( I ) = O.C 
or 20 J = 1,3 
X F ( I ) = XF(I)+R*MItJ)*X(J) 
or 21 I = 1,3 
X ( I > = XP( ! ) 

A( JA ) = A( JA ) + X ( 1 ) 

CALI PANDIH I X, I V ,RS ) 

IYT = ICOCX'RS 
I X = IV 

IF( IYT-( IYT/IP)+IP) 3C , 2 Q , 30 



- SD*C A 

SO*SA 

SO*CA 

C n*CA**?+SA **2 
<-CQ+ 1 . 0 )*SA*CA 
- SD* S A 

( -CD+ 1 . C ) *SA*CA 
CD*SA** 2 +CA **2 

1 ,NP 



2 8 


P M ( 1 , 3 ) = 


-RM( 1 , 


3 ) 




P H ( 2 , 3 ) = 


-RM2, 


3 ) 




F M ( 3 , 1 ) = 


-RM( 3, 


1 ) 




P H ( 3 , 2 ) = 


-RM 3, 


2 ) 


30 


CONTINUE 






5 C 


CONTINUE 






51 


CL NT I MUF 








DC 6C I = 


1 , NP 




6 C 


A ( I ) = A ( I ) /NT 






DC 61 I = 


I , NP 




61 


I A (I) = A ( I ) * 1 C C 
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WRITECP* 1 10) ( I A ( I ) , I = 1 ,NP) 
11C FT PM AT (IX, 3613) 

CALI DAMP (NP ,A,CEL) 

GO TO 1 

12C FORMAT < ' FLAG 1 • ) 

121 FORMAT ( • P L AG 2 • ) 

122 FORMAT ( • FLAG 3 • > 

123 FORMAT ( ' FLAG A ' ) 

END 



SLBROUT I ME 0 AMP < N P , A , DE L ) 

IMPLICIT FFAL*P < A-F » C- Z ) 

PEAL *4 A, DLL 

DIMENSION 0(10) , FINCR(IO), R(200), E ( 200 ) » I R < 200 ) 
DIMENSION F X( 200 ) , 0A(2C0), A(200) 

2C2 FCFMAT ( ' FLAG C ’ ) 

500 FORMAT ( I5,SG15.5) 

S = 0.0 

DO 10 I = 1 , NP 
S = S + n t L 

10 PX ( I ) = S 
P ( 1 ) = l.C 
P ( 2 ) = 0.0 
FINCP(l) = .00C1 
F I N'C R < 2 ) = .0001 
EFSIL = 1.0CD-C6 
DO 11 I = 1 , NP 

11 DA < I ) = A ( I ) 

CALL LEAST(2,NP,P,PX,PP,PPP,DA,FINCR,EPSIL,NGITR, 

1 R P Q » NO F , P , E ) 

PRINT 100, NO I TP , RRQ 
WP I T E ( 8 » 1 00 ) NO I TR , RRC 
IOC FORMAT ( I 5 , G 1 2 . 5 ) 

PRINT 101, ( ( P ( I) ,E(I) ), I = 1,2) 

lvPI T F(8,101) < < P< I ) ,E< I ) ) , I = 1,2) 

101 FORMAT ( IX, 1 OG 1 1.4) 

DC 20 I = 1,NP 

2C I F ( I ) = 1 0CQ*R < I ) 

Vs l PITE(R,lC2) ( I R ( I ) , I = 1,NP) 

102 FORMAT < 30 1 A ) 

RETUF N 

END 



SUBROUTINE LEAST(IR,IS,A,X,X8,XC,Z, p INCR,EPSlL,NCITR, 

1 R F 0 , NOF , P , E ) 

C 

IMPLICIT PEAL*8 (A-H,C-Z) 

C IR = NC.CF PARAMETERS, IS = N'C. OF POINTS, A IS ARRAY OF 
C PAPAMFTEFS, X IS INDEPENDENT VARIABLES CEPENOENT ) . F INCR 
C IS APFAY OF INCREMENT PARAMETERS. EPSIL IS -FRACTIONAL- 
C ERROR CPITERICN. N'CI T R IS NO. O p ITERATIONS REQUIREC (UP 
C TO 1C). PRO IS = SUM CP SQUARES OF RESIDUALS. NCF IS 
C NUMBER OF THE FUNCTION USEC IN -EON-. E IS THE ARRAY OF 
C ESTIMATED ERRORS IN THE PARAMETERS. 

C 

DIMENSION A < 10 ) , X(20C), XB(200)» XC(200 ), Z(200), FIN 
1C), 0(200,10), DT(10»200) , DEL(IO), DS(10), DP I ( 1 0 , 10 ) 
2E< 10) 

NCI TP. = 0 
DC 20 I = 1 , I S 

20 F ( I ) = 2(1) - ECN( A ,X( I ) , XR ( I ) , XC ( I ) , NCF ) 

25 IF (NOITR-F) 13C,130,A 
13C NCI TP = NCI T R + 1 
DC 220 J = 1, IR 
A( J) = A< J)+FINCR<J) 

DC 15 I = 1 , I S 

15 0(1, J) = EON ( A , X ( I ) ,Xfi(I) ,XC( I) ,NCF) 

A( J) = A ( J ) - F I N C R ( J ) 

22C CONTINUE 
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OC 3C I = 1 , I S 

CONST = FOM 4,X( I ) ,XR ( I ) ,XC(! ) ,NCF) 

DO 3C J = 1 , IR 

3C D ( I , J ) = (D( I , J )-CONST ) /F INCR ( J ) 

DO 3 5 I = It IS 
DO 35 J = 1,IR 

35 OT< J ,1) = H I , J) 

TO 36 I = I,IR 
DO 3fc J = I * I R 
DP< I » J ) = C.O 

DO 36 K = I , I S 

36 DF(IfJ) = DP(I,J)+DT(I,K)*D(K,J) 

CALL GAUSS 3 ( I R , 1 . CCE-30 ,DP , OP I , KER > 
IF(KEP-l) 120,37,120 

37 DC 4C I = 1,IS 
DC 38 J = 1 , 1 R 
DS ( J ) = C.O 

DO 38 K = 1 , IR 

38 C S ( J > = D S ( J ) + D P I (J,K)*DT(K,I) 

DO 39 L = 1 , 1 R 

3 <5 DT ( L , I ) = rS(L) 

40 CONTINUE 

DO 110 I = 1, IR 
DEL ( I ) = C.O 
DC 11C J = 1, IS 
110 DEl(I) = CEL(I )+DT( I,J)*R(J) 

DO 10 I = 1 , I R 
1C All) = A ( I ) +DE L ( I ) 

D0320 I = 1,IS 

32C P ( I ) = Z(I) - ECN( A ,X(I ) ,XB< I ) ,XC( I ) ,NCF ) 

RFC = O.C 
DO 50 I = 1,IS 
5 C RFO = RRQ+F ( I )**2 

CPES = DABSI RRC-RRP ) - EPSIL*RRP 

P t D — DC O 

IF(CRFS) 100,100,25 
101 FORMAT (2CH CONVERGENCE FAILURE) 

A PRINT 101 
GO TO 100 
1 2 C PRINT 1CC1 

1C C 1 FORMAT ( 1 6H MATRIX SINGULAR) 

ICC FIS1 = IS-IR 

DC 150 I = 1 , 1 R 

1 5 C E ( I ) = DSGRT(RRC*DPI( 1,1) ) 

RETURN 

END 



SUBFOUTINE GAUSS3 ( N , EP , A , X ,KER ) 

IMPLICIT RFAL*8 (A-H, C-Z) 
DIMENSION 4(10,10), X(1C,10) 

DC 1 1=1, N 
DO 1 J= 1 , N 

1 X(I,J)=0.0 
DC 2 K= 1 , N 

2 X ( K , K ) = 1.0 

10 DO 34 L = 1 , N 
KP = 0 

Z = 0.0 

DC 12 K=L , N 

I F ( Z-DAB S ( A ( K , L ) )) 11,12,12 

11 Z=DAPS( A ( K ,L ) ) 

K p — K 

12 CONTINUE 
IF(L-KP)13,20,2C 

13 DC 14 J = L , N 
Z = A ( L , J ) 

A(L,J)=A(KF,J) 

14 A ( KP , J ) = Z 
DO 15 J=1,N 
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7 = X ( L » J ) 

X(L,J)=X(KP,J) 

15 X ( KP , J ) = 7 

2C IF(PABS(A(L,L) ) -E P ) 50 , 5 C , 30 

30 I F ( L-N ) 31 » 34 ,34 

31 L P 1 = L + 1 

DO 36 K=LP1,N 
I F ( A < K , L » ) 32 , 3 6 , 32 

32 RATIO=A(K,L)/A(L,L) 

PC 33 J=LP1,N 

33 A(K,J)=A(K ,J)-RATIC*A(L»J) 

PC 35 J = 1 » N 

35 X(K, J)=X(K, J)-RATIC*X(L,J) 

36 CONTI NUE 

34 CC NT I NUF 

40 DO 43 1 = 1, N 
I I =N + 1- I 

DO 43 J = 1 , N 
S = 0.0 

I F ( I I-N )41 ,43,43 

41 IIP1=II+1 

DO 42 K = I I P 1 , N 

42 S = S + A ( l I ,K )*X< K, J ) 

43 X< II , J » = ( X< I I , J ) - S ) /A( I I , I I ) 

KE R= 1 

RETURN 
50 KEP=2 
RETURN 
END 

REAL FUNCTION E CN 4 8 ( A , X , X B , XC , N CF ) 
IMPLICIT REAL*« (A-H,C-Z) 

DIMENSION A ( 10) 

EON = PEXP(-A(2)*X)*CCGS(A( 1 )*X) 

PFTURN 

END 
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